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The roles of isospin asymmetry in nuclei and neutron stars are investigated using a range 
of potential and field-theoretical models of nucleonic matter. The parameters of these models 
are fixed by fitting the properties of homogeneous bulk matter and closed-shell nuclei. We 
discuss and unravel the causes of correlations among the neutron skin thickness in heavy 
nuclei, the pressure of beta-equilibrated matter at a density of 0.1 fm -3 , the derivative of 
the nuclear symmetry energy at the same density and the radii of moderate mass neutron 
stars. Constraints on the symmetry properties of nuclear matter from the binding energies 
of nuclei are examined. The extent to which forthcoming neutron skin measurements will 
further delimit the symmetry properties is investigated. The impact of symmetry energy 
constraints for the mass and moment of inertia contained within neutron star crusts and the 
threshold density for the nucleon direct Urea process, all of which arc potentially measurable, 
is explored. We also comment on the minimum neutron star radius, assuming that only 
nucleonic matter exists within the star. 
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I. INTRODUCTION 

Strongly interacting matter in which more neutrons than protons exist is encountered in both 
heavy nuclei and neutron stars. In stable nuclei the net asymmetry / = (N — Z)/(N + Z) ranges 
up to about 0.24, but the neutron-to-proton asymmetry, 5 = (n n — n p )/n where n n and n p are 
the number densities of neutrons and protons and n = n n + n p , approaches unity in the nuclear 
surface. In the future, rare-isotope accelerator experiments will extend the range of / to values 
well in excess of 0.2. In contrast, / could be as large as 0.95 in the interiors of neutron stars. 
The physical properties of nuclei, such as their masses, neutron and proton density distributions 
(including their mean radii), collective excitations, fission properties, matter and momentum flows 
in high energy heavy-ion collisions, etc. all depend on the isospin structure of the strong interactions 
between nucleons (i.e., nn and pp interactions versus np interactions). The energetics associated 
with the n — p asymmetry can be characterized by the so-called symmetry energy, E sym , which is 
the leading coefficient of an expansion of the total energy with respect to asymmetry: E(n, 5) ~ 
E Q (n)+E sym S 2 ---. The energy fi = [i n — u p = AE sym 5, where \i n and \i p are the neutron and proton 
chemical potentials, respectively, is crucial in determining reaction rates involving electrons and 
neutrinos, particle abundances, etc., in astrophysical contexts such as supernova dynamics, proto- 
neutron star evolution, the r— process, the long-term cooling of neutron stars, and the structure of 
cold-catalyzed neutron stars (i.e. their masses, radii and crustal extent), etc. The pervasive role of 
the isospin dependence of strong interactions in nuclear processes in the laboratory and the cosmos 
is sketched in Fig. 1. In this work some of these connections will be discussed. 
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FIG. 1: The multifaceted influence of the nuclear symmetry energy. 

Recently, several empirical relationships have been discovered that underscore the role of isospin 
interactions in nuclei and neutron stars. These include correlations between: 

1. The neutron star radius R and the pressure P of neutron- star matter. Lattimer and Prakash 
0, 0| found that the quantity i?P _1//4 is approximately constant, for a given neutron star 
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mass, for a wide variety of equations of state when the pressure P of beta-equilibrated 
neutron-star matter is evaluated at a density in the range no to 2no, where no denotes 
equilibrium nuclear matter density. Since the pressure of nearly pure neutron matter (a 
good approximation to neutron star matter) near no is approximately given by n 2 dE sym /dn, 
the density dependence of the symmetry energy just above no will be a critical factor in 
determining the neutron star radius. 

2. The neutron skin thickness in nuclei and the pressure of pure neutron matter at sub-nuclear 
density: Typel and Brown 0, have noted that model calculations of the difference between 
neutron and proton rms radii 5R = (r 2 ) 1 ^ 2 — (rp 1 } 1 ^ 2 are linearly correlated with the pressure 
of pure neutron matter at a density below no characteristic of the mean density in the nuclear 
surface (e.g., 0.1 fm -3 ). The density dependence of the symmetry energy controls 5R (we 
will call this the neutron skin thickness) in a heavy, neutron-rich nucleus. Explicitly, 5R 
is proportional to a specific average of [1 — E sym (no) / E sym (n)] in the nuclear surface, see 



jporti 

Refs. 0, and Eq. (|TTf2l below. 

3. The neutron skin thickness in nuclei and the neutron star radius: Horowitz and Piekarewicz 
01 have pointed out that models that yield smaller neutron skins in heavy nuclei tend to 
yield smaller neutron star radii. These authors, along with others @, 0], have also pointed 
out the need for an accurate measurement of the neutron skin. 

Unlike proton distributions, neutron distributions in nuclei have remained uncertain to this date. 
Recent studies of neutron densities from a global analysis of medium-energy proton scattering on 
208 Pb indicate that 0.07 < 5R < 0.16 fm [ic| . Related information is also available from an 
analysis of antiprotonic atom data that gives 5R = 0.15 ± 0.02 fm jlH ]. in the latter work, nucleon 
density distributions are parameterized by Fermi functions and it is found that the half-density 
radii for neutrons and protons in heavy nuclei are the same, but the diffuseness parameter for the 
neutrons is larger than that for the protons. Skin thicknesses as large as 0.2 fm were obtained 
in earlier analyses; see the discussion of Karataglidis et al. ^3). Since these studies involve 
strongly interacting probes, even to this date the value of SR for a nucleus such as 208 Pb is not 
accurately known. This situation should improve as it is expected that the neutron rms radius will 
be determined to about 1% accuracy by measuring the parity- violating electroweak asymmetry in 
the elastic scattering of polarized electrons from 208 Pb an experiment planned at the Jefferson 
Laboratory 

A common theme in the evaluation of the neutron skin thickness and the pressure of neutron 
star matter below and above the saturation density no is the density dependence of the isospin 
asymmetric part of the nuclear interaction. To highlight this dependence, we examine theoretical 
predictions of 5R and the neutron star radius for an extensive set of potential and field-theoretical 
models. For the former we employ Skyrme-like potential model interactions, including a param- 
eterization [lii ] of the microscopic calculations of Akmal and Pandharipande ■ For the latter 
field-theoretical models we use Walecka-type models (lj| in which nucleons interact through the 
exchange of a, u> and p mesons augmented with additional mesonic couplings (see, for example, 
Horowitz and Piekarewicz 0]) to better describe the properties of nuclei and explore the properties 
of neutron stars. In order to provide baseline calculations, our investigations will be restricted to 
the case in which the high density phase of neutron stars contains nucleons (and enough electrons 
and muons to neutralize the matter) only. 

Our goals in this work are to identify the causes of the correlations mentioned above through a 
semi-analytic analysis of detailed calculations and to uncover other possible correlations. Toward 
these goals, we investigate an assortment of potential and field-theoretical models that cover a 
range of symmetry properties and incompressibilites, but which are constrained by the empirical 
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properties of isospin symmetric and asymmetric nuclear matter, and the binding energies and radii 
of closed shell nuclei. The properties of nuclei are investigated through Hartree-Fock-Bogoliubov 
calculations for potential models and Hartree calculations for field theoretical models. In both 
categories, our calculations include models that match the Akmal and Pandharipande equation 
of state (EOS) as results for nuclei with this model are not available in the literature. Additional 
insights are provided through analytical and numerical analyses of isospin asymmetric semi-infinite 
matter in the potential and field theoretical approaches. The results of these calculations are 
utilized in understanding several correlations in conjunction with inferences from nuclear masses. 
The neutron skin-pressure and neutron star radius-pressure correlations are reexamined to provide 
improved relations. Our calculations also provide an estimate of the smallest neutron star radii 
that ensue from these models. 

The presentation is organized as follows. In Sec. II, the Hamiltonian (Lagrangian) densities of 
the potential (field-theoretical) models used in this work are described. The variational approach to 
determine the bulk and surface properties of isospin asymmetric semi-infinite nucleonic matter in 
these models is developed in Sec. III. This section also contains a discussion of how these properties 
enter in the liquid droplet model description of nuclei. Our results and discussion of the various 
correlations, and their origins are presented in Sec. IV. Section V contains a brief description of 
other related correlations. Discussion and conclusions are contained in Sec. VI. Appendix A lists 
the properties of the models used in this work. The coupling strengths of some newly-constructed 
models are given in Appendix B. 



II. EFFECTIVE THEORIES 



In both potential and field-theoretical model descriptions of many-body systems, the Hohenberg- 
Kohn-Sham theorem [lti H?) ] allows the total Hamiltonian density to be expressed in terms of the 
local particle number and kinetic energy densities of the various species. This simplification is of 
great utility in the study of both heavy nuclei and infinite matter since microscopic calculations, 
such as variational/Monte Carlo studies starting from a Hamiltonian constructed on the basis of 
nucleon-nucleon and many-nucleon interactions, are computer intensive and are only beginning 
to be undertaken for heavy nuclei (for a review of such studies on light nuclei, see Ref. [la]). In 
practice, the strengths of the various interactions in the density functional approach are determined 
by appealing to available data. In the following, we consider both potential and field theoretical 
approaches so that distinguishing traits between these approaches can be identified. 



Potential Models 



A commonly employed non-relativistic effective Hamiltonian density stems from Skyrme's work 
in which finite-range forces between nucleons were approximated by effective zero-range forces. 
For low relative momenta (i.e., with s— and p— wave interactions only), the resulting Hamiltonian 
density takes the form [2(3] 

(1) 



H = H B + \ [Qnn(^n n ) 2 + 2Q np Vn n ■ Vn p + Q pp (Vn p ) 2 \ +H C + Hj , 
where the term associated with spatially homogeneous bulk matter is of the form 



Hb 



Tr 



2m, 



T n + 7, T P +n{T n + T p ) 



+ {j n n n + T p n p ) 



t(] +X2 )-tG +Xi 



2 

t 3 



2 



+ 



12 



1 + 



X3 



n 



n 



+ x 3 



2 , 2 

n„ + n p 



2 i 2 

n n + n p 



n 



(2) 



Above, m n and m p are the masses of the neutron and the proton (these will often be taken to be 
equal and then denoted by M), n n and n p are the number densities of neutrons and protons, and 
to, t±, t2, t 3 , xq, xi, X2, x 3 , and e are parameters that give the strengths of the various potential 
interactions. We recall that total baryon density is n = n n + n p . The first two terms represent the 
kinetic energy densities for the neutron and the proton respectively, the sixth term is the density- 
dependent form of the zero-range multi-body force, and the remaining terms constitute the s— and 
p— wave parts of the two-body interaction. 

The coefficients Q associated with the spatially varying part of the Hamiltonian density in 
Eq. (J2) are given by 
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Note that for a Skyrme-type force, the Q's are constants. However, in general, it is possible for the 



Q's to be density-dependent functions. Also note that, for Skyrme-type forces, Q r 
this does not have to be true in general. Generally, we have Q np = Q pn however. 
The term Tic hi Eq. arising from the Coulomb interaction is 



n c {f\ 



e 2 n p (r) 



3e 2 



1/3 



n 



,(r) 4/3 , 



Q pp , but 



(5) 



where e denotes the proton's electric charge and the second term is an exchange correction |2 
The term Tij arising from the spin orbit interaction is 

_Wo 
2 
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(n n V ■ J n + n p V ■ Jp + raV 
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(6) 



ip-a x vipi, similarly for protons and J 



where the neutron spin-orbit density J n 

Jn ~\~ Jp- 

The use of Skyrme's effective interactions has successfull y re produced many of the global prop- 
erties of nuclei including binding energies, charge radii, etc. |22j. The properties of bulk symmetric 
nuclear matter at the saturation density no, can be expressed directly in terms of no, and the 
model parameters. Explicitly, the energy per nucleon E/A, the effective mass M*, the pressure 
per particle P/n and the compressibility K are given by 




1 + 0no) + —n + —^o 
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The symmetry energy of bulk matter, E sym , is 



n 2 d 2 E/A 



2 da 2 
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(12) 



where the isospin asymmetry density is denoted by a = n n — n p . The volume symmetry energy 
of equilibrium nuclear matter E syrn (rio) will be denoted S v . We note that the symmetry energy 
asymptotically decreases at large densities if e > 2/3 or t\X\ > t2(4 + 5x2) /3. (This poor high- 
density behavior can be remedied H^|.) For models in which the symmetry energy monotonically 
increases with density, the negative sign of the multi-body term involving t3 generally results in 
neutron star radii which are smaller than in field-theoretical models (see Sec. IIV C 3|) . 



B. Field Theoretical Models 

The structure of nuclei and the properties of high-density nucleonic matter have also been 
studied utilizing the mean field approximation to a Walecka-type Lagrangian that couples nucleons 
to scalar (a), vector (u^), and vector-isovector (p At ) mesons and photons (A^) [3]. The Lagrangian, 
supplemented by non- linear interactions, is 



C 



-V(a) - If^r + Iml^u,* ~ ■ + WJ» ■ P, - \F IIV F^ 



4- 1 fiiy J 



24 



where f are the SU(2) isospin matrices, f^, v 
and the scalar meson potential 



V{a) = \mW + ^{g„af + ±{ g<J af 



(13) 

d^Ay Q v A p 

(14) 



In Sec. IV, we discuss the choice of the couplings g a , g u , g p , k (dimensionful) and A, as well as £ 
and £ that allow the high-density behavior of the equation of state to be varied. The first five of 
these are constrained to reproduce the properties of equilibrium nuclear matter and the empirical 
symmetry energy. Horowitz and Piekarewicz 0] have extended the non-linear Walecka model by 
adding a 2 p^ • p p and uj^u^p^ • p^ terms to the Lagrangian so as to modify the density dependence 
of the symmetry energy at supranuclear densities (in the absence of this coupling the dependence 
is linear at large densities). In order to provide additional freedom in varying the symmetry energy 
we have adopted a general function /, namely 



(15) 



Utilizing the mean field approximation with the Lagrangian 1)13(1 , the binding energy per particle 
in symmetric matter at equilibrium is given by 



B = M 
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where the Dirac effective mass is M* = M — g a ao, and the Fermi momentum and energy of 
symmetric nuclear matter are, respectively, 



3vr 2 n 
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(17) 



The equilibrium pressure per particle is 
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and the compressibility at equilibrium is given by 
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The symmetry energy of bulk matter is 



k 2 

E S ym(n) = 77~~" + 
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QE* F 8mM + 2/K,^ 



(19) 



(20) 



In the above equations, the subscript "0" denotes mean field values of a and lj. For the case / = 0, 
the symmetry energy in Eq. (|2Uj) varies linearly with the density at large densities. The presence 
of the function / permits variations in this density dependence. 



C. The EOS of Akmal-Pandharipande-Ravenhall 



The results of the variational microscopic calculations of Akmal and Pandharipande |lj], in 
which many-body and special relativistic corrections are progressively incorporated into prior mod- 
els, have been parameterized by Akmal, Pandharipande and Ravenhall (hereafter APR) The 
bulk effective Hamiltonian density is 

Hp, APR = (j^ + (P3 + {l-x) P5 ) ne- p4n ^j r n + + (p 3 + xp 5 ) ne" P4n ^J t p 

+ 9l (l - (1 - 2x) 2 ) +92 (l-2x) 2 , (21) 



where r n and t p are the kinetic energy densities of the neutrons and protons, respectively and 
x = rip/n. Contributions from potential interactions are contained in the functions g\ and g2 that 
have two different forms corresponding to the low density phase (LDP) and the high density phase 
(HDP) that contains a neutral pion condensed phase. Explicitly, 

giLDP = -n 2 (pi + p 2 n + p 6 n 2 + (pio + pure) e~ p 9 n2 ) , (22) 
g2LDP = -n 2 (p 12 /n + p 7 +p$n+p 13 e~ p v n2 ^ , (23) 
9ihdp = 9ildp - n 2 (p 17 (re - p 19 ) + p 2 i (n - Pl9 ) 2 ^"(""Pm)) , (24) 
g2HDP = g2LDP ~ n 2 (p ls (re - P20) + Pu (re - P20) 2 e pl ^ n ~ P20 ^ . (25) 
The numerical values of the various parameters pi for the recommended "A18+fo-|-UIX*" EOS, 



9 



TABLE I: Coupling strengths for the Hamiltonian density of APR in Eq. i|21|) . The dimensions are such 
that the Hamiltonian density is in MeV fm -3 . 



Pi 


P2 


P3 


Pi 


P5 


P6 


P7 


Pa 


P9 


Pio 


Pn 


337.2 


-382 


89.8 


0.457 


-59.0 


-19.1 


214.6 


-384 


6.4 


69 


-33 


Pl2 


Pl3 


Pu 


Pl5 


Pie 


Pl7 


Pis 


Pl9 


P20 


P21 




0.35 








287 


-1.54 


175.0 


-1.45 


0.32 


0.195 








that incorporates three-nucleon interactions and relativistic boost corrections, are collected in 
Table HI For this EOS, n = 0.16 fm~ 3 , E/A = -16 MeV, M*/M = 0.7, K = 266 MeV, and 

VW = 32 - 6 MeV - 

For later use the Q's for the EOS of APR, calculated using the procedure detailed in Ref. |2a |. 

are 

- P4(P3 ~ 2p5)(ra n + 2n p )] , 

- P4(P3 - 2p5)(ra p + 2n n )] , 

4p 5 ) - 3p 4 (p 3 - 2p 5 )(n n + n p )] . (26) 

Note that Q nn / Q pp except in the case that n n = n p . 

A Hartree-Fock calculation of finite nuclei using the bulk Hamiltonian density in Eq. (|21j) aug- 
mented by gradient, Coulomb and spin-orbit terms is more complicated than that using the Skyrme- 
like Hamiltonian in Eq. (J2J) because of terms that vary exponentially with density in Eq. (|21|) . 
Therefore, we have fit the APR equation of state for bulk nuclear and neutron matter up to 3/2 
nuclear matter density to that of a Skyrme-like Hamiltonian. For a realistic description of nuclei, 
the density dependence of the effective masses is important, so we have constrained the fit so that 
the effective masses match as closely as possible those of the APR equation of state. We have 
also adjusted the coefficient of the spin-orbit interaction so as to closely match the charge radii 
and binding energies of 208 Pb, 90 Zr, and 40 Ca. This procedure leads to the parameter set given in 
Table HT1 and is referred to as NRAPR below. Figure |2] shows that in matter a good fit to the exact 
APR binding energy /particle is obtained. A fit was also made using the field-theoretical approach, 
referred to as RAPR, and the parameters are tabulated in Table III. As shown in Fig. [2 the fit is 
even a little better than in the non-relativistic case. 



TABLE II: Parameters for the Skyrme Hamiltonian obtained from a fit to the EOS given by the APR 
Hamiltonian, Eq. I|21(l ; the dimensions are such that H is in MeV fm~ 3 . This model is referred to as 
NRAPR in the text. 



to 
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t 3 


Xq 


-2719.7 
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15042 
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Xl 


X2 
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-0.047986 


0.027170 


0.13611 


0.14416 


41.958 



It is difficult to utilize these low-density fits to APR at densities of relevance to neutron stars, 
largely because of the phase transition that is present in APR. There is not enough freedom in either 
the present potential model approach or the field-theoretical approach to simultaneously match the 
equation of state at low densities and satisfy the constraint on the maximum mass neutron star. We 
have also attempted fitting the low-density and high-density phases of APR separately. Although 
this is somewhat difficult to accomplish with the field-theoretical model above, it is easily done 



Qnn — 




-p 4 n 


-6f> 5 


Qpp — 


¥ 


-p 4 n 


-6f> 5 


Qpn 




-p 4 n 


[4(p 3 ■ 



10 




FIG. 2: Binding energy per particle versus density in nuclear matter and neutron matter for the EOS's of 
APR (solid line), NRAPR (dashed line), and RAPR (dotted line). 



TABLE III: Coupling strengths for the Lagrangian in Eq. 1|13|) that fit the APR Hamiltonian density in 
Eq. I|21|) up to 3/2 nuclear matter density; the dimensions of a, and bi are such that £ is in MeV 4 . This 
model is referred to as RAPR in the text. 
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9p 
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494.28 MeV 


7.8144 


9.2629 


10.288 


4.3789 MeV 


0.078205 


c 
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cti 


do 


a<3 


(24 


0.065381 


0.30524 


0.18566 


0.73588 


7.5995 xl0~ 3 


3.4620 x!0~ 4 


a 5 


a 6 


h 


b 2 


b 3 




6.7880 xl0~ 6 


5.4824 xl0~ 8 


0.039169 


3.1323 xl0~ 9 


9.8907 xlO" 10 





with the potential model. However, this patchwork approach falls outside our goal of studying 
correlations that are generic to equations of state without softening phase transitions. In what 
follows, therefore, results for nuclei are calculated using the NRAPR and RAPR models, whereas 
those for neutron stars are from calculations using the original EOS of APR. 



III. SEMI-INFINITE NUCLEONIC MATTER 



The isospin dependence of the nuclear force has its most important ramifications in the bulk 
and surface energies of nuclei. In finite nuclei, these effects are coupled, but certain properties 
like the neutron skin thickness are affected more by the surface effects. The role of the surface is 
most easily ascertained by an investigation of semi-infinite matter 2f], 2^. In this case, matter is 
assumed to vary only along one axis (the z-axis) and is assumed to be uniform in the two remaining 
spatial directions. Matter at the z — > — oo limit is saturated matter at a specified proton fraction, 
whereas the z — > +oo limit is the vacuum (see Fig. Matter in these two limits are in chemical 
and pressure equilibrium with each other. Thus, matter in the z — > — oo limit is at the saturation 
density for the reference proton fraction, with vanishing pressure. 
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FIG. 3: Schematic diagram of the surface regions of semi-infinite matter indicating the relation between the 
neutron skin thickness, SR, and the corresponding quantity t = R n — R p determined from the squared-off 
neutron and proton density distributions. 



We note that in sufficiently neutron-rich matter, such as neutron star matter, the (non- 
relativistic) chemical potential for neutrons at the saturation density is positive. In this situation, 
which occurs when the reference proton fraction is less than about 0.35, the pressure at the satu- 
ration density will also be positive and there must be a finite density of neutrons in the z — > +oo 
region with the same pressure and neutron chemical potential. 

We will study how the surface energy and the neutron skin thickness vary in the semi-infinite 
approximation due to modifications of nuclear model parameters (including those that delineate 
isospin effects). These features can be connected to experimental properties of laboratory nuclei. 



A. Variational Approach 

The density profiles of semi-infinite matter (schematically illustrated in Fig. ^ are obtained 
by minimizing the Hamiltonian (energy) density subject to the constraints of baryon number and 
charge conservation. The baryon and isospin asymmetry densities, 

n = n n + n p , a = n n - n p , (27) 

respectively (and fields and functions of them), are z-dependent even though we do not explicitly 
indicate them as such. The neutron excess is defined by 5 = a/n. We will denote quantities at 
z = — oo by the use of the subscript "L" . Primes will generally indicate derivatives with respect to 
the coordinate z. Pressure equilibrium between the dense nuclear phase and the vacuum requires 
that the pressure Pl = 0. Therefore, the density ni is the saturation density for matter with a 
neutron excess 5l- The value of n^, the surface profiles n n and n p , and the surface energy (to be 
defined below), are all determined for a given nuclear force model by the single quantity 5l- 
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We denote the saturation density of symmetric matter by no = til(5l = 0). The symmetry 
energy function E sym (n) is generally a function of density and temperature T, but we consider 
only T = matter in this paper. For isospin symmetric matter E sym (no) = S v is the standard 
volume symmetry energy parameter whose value lies in the range 25-35 MeV. We also define the 
squared-off neutron and proton radii (see Fig- El) 7 Rn an d R P , by 



OO oo 

n n L{Rn - L) = J n n dz ; n pL (R p - L) = J n p 



dz . 



(28) 



where L — > oo. Although these integrals formally diverge, the quantities R n and R p are finite. 
The variation of the Hamiltonian density 



oo 

5 J \H — un — fA a a] dz 



(29) 



results in the two Euler-Lagrange equations 



dn 

dn 

an 

da 



d dn 

^ dz dn' 

d dn 

^ a dz da' 



(30) 
(31) 



where 2u = [i n + Hp and 2[i a = fi n — fi p are the Lagrange multipliers (equivalently, chemical 
potentials) associated with the constraints of baryon number and charge conservation, respectively. 
These chemical potentials are related to the individual neutron and proton chemical potentials 
through fin + [x a a = ^ ri n n + fi p n p . 

For Hamiltonians of the form of Eq. (^Q), without 7^c and Wj, we find 



dn B 

dn 



1 1 



&Hb 
da 



He 



n" , 
T (Qr 

1 

+ 77 



a 

2Qnp Qpp) H {.Qr 



Q 



dQ 



Tip 



+2 



dn 
( dQ, 
\ da 



(V) 2 + («0 2 ) 



(dQ 



n 2 ^Qnp _|_ 

da 



V dn 
dQpp\ __/_/ 



ppj 
n . dQ pp 



da 



n'a + 2 



dn 
( dQ r , 



a" 

T (Qr 



n 



2Qnp Qpp) , \Q 



1 



dQ 



Tip 



da 

(dQnn 
dn 



l\2 i / l\2\ , ( QQr 



n 



+ W) + 



V da 



2 dQnp + dQpp | n , a 



dn 



dn 



\ da 

Qpp) 

dQpp 
da 

(dQnn 
dn 



{n'f - 

dQ pp 
da 



(a 



>\2 



a 



+ 



{{n'f - (a') 2 ) 
dQ pp 



dn 



(32) 



(33) 



Both sides of these equations vanish for z — > ±oo. 

The surface energy is the difference between the total energy and the bulk energy of an equivalent 
number of neutrons and protons at the saturation density. Equivalently, the surface thermody- 
namical potential per unit area is expressed [3] by the same integral that is minimized in Eq. 
CU): 



a = J [n — un — u a a] dz . 



(34) 
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The quantity a is also known as the surface tension. The integrand vanishes at z — > ±00, expressing 
the fact that the pressure vanishes there. Expanding <j{8l) in terms of the neutron excess 5l 



a(S L ) = a{5 L = 0) + ^ 

L 



5i + ... = a - a s 5 2 L + ... (35) 

s L =o 



allows us to identify a$ as the symmetry parameter of the surface tension. We note that the surface 
energy (which is the product of surface tension and a characteristic surface area) , as distinct from 
the surface thermodynamic potential, has a surface symmetry parameter identical to this, but of 
opposite sign |2^ |. 

The neutron skin thickness is defined by 

t = R n -R p = 7 ( - 2l) dz . (36) 
J \ n nL n pL I 

—00 x ' 

Using the relations 

n-nL = \n L (1 + 5 L ) ; n pL = \n L (1 - S L ) (37) 

the neutron skin thickness can be expressed in the equivalent form 

00 

* = :-T7^T / n(l-l)dz. (38) 



n L (l-5 2 L ) J \S L 

—00 

Were the neutron and proton density distributions describable by Fermi functions defined by half- 
density radii z UtP and diffuseness parameters a n ^ p . 

n in,p)L (39) 



the skin thickness would be t = z n — z p , to lowest order in the diffuseness parameters a np . We 
note that the definition of skin thickness might be problematic in finite nuclei where there are 
often oscillations of the density profiles near the center. In this case, the skin thickness is better 
expressed as the difference between the neutron and proton rms radii (r^p) 1 ^ 2 . To lowest order in 

diffuseness corrections, one finds that 5R = (r^) 1 ^ 2 — (r 2 ) 1 ^ 2 = i/| t. 



1. Potential Models 

The Euler-Lagrange equations (|3*2|) and ()33|) . multiplied respectively by n' and a' , then added 
and integrated, yield 

in') 2 (a 1 ) 2 n'a! 

{Qnn ^Qnp ~i~ Qpp) q h {Qnn ^Qnp ~i~ Qpp) „ H {Qnn Qpp) "j — 7~(-B ~ 

fxn - fi a a . (40) 

o o 4 

This equation holds even in the case in which the Q's are density-dependent. The right-hand 
side of this equation is the bulk thermodynamic potential density, Qb/V = —Pb, where Pb = 
n 2 [d{H.B /n) /dn] is the bulk pressure. Note from Eq. ()40j) that the bulk and gradient terms 
contribute equally to the surface tension in Eq. (|34|) . which can therefore be written as 

00 

<j(5l) = 2 J [Hb — fJ-n — fj, a a]dz . (41) 

— oo 
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To lowest order the bulk Hamiltonian should be quadratic in the neutron excess 5: 

H B = n B (n,5 = 0)+E sym (n)n5 2 . (42) 
Since the bulk pressure vanishes for z — ► — oo, the chemical potentials, to second order in 6l, are 



H n = -B + S V 5 L (2 - 5 L ) , 
fi p = -B- S V 5 L (2 + 5 L ) , 



(43) 
(44) 



where B = -Hb(wo, S = 0)/no is the symmetric matter bulk binding energy per particle. Therefore 
the surface tension becomes 



H B (n, 6 = 0) + n[B + 5 z E sym (n) + 5 Z L S V - 255 L S V 



dz . 



(45) 



We are interested in evaluating this in the case that 5l is small so that we can take Qnn — Qpp] this 
is exactly true for standard Skyrme Hamiltonians and approximately so for the APR Hamiltonian. 
It is also reasonable in this case to use Eq. (001) with (a') 2 set to zero which results in a quadrature 
formula for the surface density profile: 



Qnn + Qnp dn = 2^JHb - fin - fi a a dz . 
Utilizing this result, Eq. (|45|) can now be written as 



(46) 



" L 



ct(5 l ) = y/Q nn + Qnp I [n B (n, 5 = 0) + n (b + 5 2 E sym (n) + 5 2 L S V - 255 L S^] 1/2 dn . (47) 



Again setting Q nn = Q pp , neglecting derivatives of a and further assuming that the dependence of 
the Q's on a can be neglected, the right hand side of Eq. is zero so that 



AE sym (n)6 ~ n n - u p . 
Using Eqs. (|4^|) and (|4"4"|) this gives the useful result 

S S v 



5l E sym (n) 



(48) 



(49) 



(Note that the use of Eq. (|48|). which is strictly valid to order 5 in bulk matter, requires that the 
variation of 5 in the surface tracks its behavior in the bulk.) Substitution of these relations in Eq. 
(|4*7|) yields the quadrature 



C"(5l) = \J Qnn + Qnp j 



H B (n,5 = 0) + nB - nS v S 2 L 



S v 



Esym ip>) 



1/2 



dn . 



(50) 



Then the surface tension of symmetric semi-infinite matter is 



oo = ct(Sl = 0) = Q nn + Q n p J [Tt B (n, 5 = 0) + nB] 1/2 dn . 



(51) 
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and expansion of the argument of the integral in Eq. (|50[) allows us to write the surface symmetry 
parameter in the convenient form 



I'L 



Q 



up 



n 



S v 



[H B (n,5 = 0) + nB]- 1/2 dn. 



(52) 



Our results will be given in terms of the surface energy E s = 4-7rrQ<7o and the surface symmetry 
energy S s = 4irr 2 a$, where (47rrQ/3)(0.16 fm~ 3 ) = 1. It is often useful to eliminate the Q's by 
forming the ratio of Eqs. (|52j) and H51|) . This yields 



Ss 

S v 



J u [S v /E S y m \ 
o 



n 



1] [H B (n,5 = 0) + nB] 



-1/2 



dn 



"L 

.! 





J [H B {n,5 = 0) + nB] 1/2 dn 



(53) 



This equation and Eqs. ()51j) and (|52*)) are central to the results to be presented in Sec. IV below. 
Finally, the profile equation (|4Ujl and the relation for 5, Eq. (|49|) . can be employed in the expression 
l|3*5|) to yield the relation 



25 L 



<?5 



(54) 



n L (1 - 5\) S v 

between the skin thickness and the surface tension that was conceived in the droplet model [2^ . 
3ol . 3ll | and later shown 0] to be quite general. 

The analytical description of the nuclear surface based on the leptodermous expansion was 
pioneered by Myers and Swiatecki |2{| and developed further in several works |R Eil. liol EH. I32T] . 
Our analysis here offers a convenient expression for S s /S v (Eq. (|53|)) in terms of the energy density 
functional for isospin symmetric infinite nucleonic matter. 



2. Field Theoretical Models 

For spatially non-uniform semi-infinite isospin asymmetric matter, the Euler-Lagrange equations 
corresponding to the Lagrangian ()13j) yield the meson and baryon field equations 





= rnlao 




= ml^o 


Pi 


= rn 2 p Po 










= 



.-g a a + -g a a -g pPo — 

C 4 3 1 2 2 ®f 



(55) 

(56) 

(57) 
(58) 
(59) 

where primes indicate derivatives with respect to z, the subscript "0" indicates that the mean field 
approximation has been used and po refers to the mean field of the neutral rho meson. We shall 
adopt the Thomas-Fermi approximation since it provides a reasonably accurate simplification of 
the Hartree approach [2^]. Here the meson fields are assumed to vary sufficiently slowly that the 
nucleons can be regarded as moving in locally constant fields at every point in space. Thus the 
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Fermi momenta and the effective mass depend upon z, i.e., kp n (z), kF p (z) and M*(z). The local 
neutron, proton and scalar densities at a given z are then 

_ Mk %l 

nn - 3vr2 ' Hp ~ 3tt 2 ' 

n - 2 f — 2 / M * 

Us ~ J (2tt) 3 V/k 2 + M* 2 + 7 (2tt) 3 VA: 2 + M* 2 ' ^ ' 

respectively. The energy eigenvalues of the Dirac equations for the neutron and proton fields are 



e n (k) = Vk 2 + M* 2 + - |5pPo 5 e p (A;) = Vk 2 + M* 2 + ^u; + ± 9p p , (61) 
respectively. The local Thomas-Fermi Hamiltonian density is then 

1 f 12 , „2 2\ 1 / ,/2 , _2, ,2\ 1 ( J2 , ™2„2\ , K / _ \3 . ^ /„ _ \A 



W(z) = - [a' 2 + m l a o l ) - - (u' 2 + m 2 ^ 2 ) - - [p 12 + m 2 pP 2 ) + - (g a a f + — (g„<T 



k Fn k F P 

C / \ 4. C / \ 4 22^ f d k f d k 

■24(5^^0) -^(sppo) -5 P Po/ + 2 y (^)3 en + 2 y ~(2^f ep 



Note that, in contrast to the nonrelativistic case, the energy density is not simply expressible in 
terms of just the baryon densities and their gradients. The ground state of the system is obtained 
by minimizing Tt(z) — fi n n n — fjL p n p with respect to the local Fermi wavenumbers kp n (z) and kF p (z). 
The Lagrange multipliers (chemical potentials) that fix the number of neutrons and protons are 
then 

Hn = e n {k Fn ) ; fi p = e p (k Fp ). (63) 

The algebraic nature of these two coupled equations relating the neutron and proton density profiles 
to the meson fields stems from the Thomas-Fermi approximation. 

In order to obtain the density profiles n n (z) and n p (z), Eqs. (|55j) - (|57|) and (|63j) are supplemented 
by the boundary conditions 



pq{-L) = <tql, u (-L)=u} 0L , Po(-L) = Pol 
a (L) = 0, cu (L)=0, Po(L)=0, 



L — > 00 



where the subscript "OL" indicates values of the fields in equilibrium at the specified neutron excess 
5l- (Note that the meson fields are finite, but exponentially small, in regions of space where the 
baryon density vanishes as prescribed by the Thomas-Fermi approximation.) These conditions 
guarantee that the derivatives of the fields vanish at the boundaries. In practice, we impose these 
conditions at points (— L,L) that are sufficiently far from the surface to ensure that the surface 
energy is unchanged by small variations in L. 

The numerical solution of the differential equations requires some care because of their behavior 
at the boundaries. We find it best to use a relaxation method. Starting with small values for the 
field derivatives, we integrate from —L outwards. The integration has to be stopped when one of 
the fields becomes negative, say at L\ > —L. This solution can then be used as an initial guess 
for relaxation in the region (— £2,-^2), where Li = (L\ + L)/2. After convergence is obtained, L2 
can be increased by a small amount and the relaxation applied again until the field derivatives are 
sufficiently small to ensure that the surface tension remains stationary. 

The determination of the density profiles enables the calculation of the surface tension. In 
the non-relativistic case Eq. (|4U|) showed that the bulk and gradient terms in the Hamiltonian 
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density gave equal contributions to the surface tension. A similar result can be obtained here by 
multiplying Eqs. (|55|) - (|57|) by the derivative of the appropriate field with respect to z, combining 
the equations and integrating. This gives 

H B - Am - fJ- a a = \ {o'$ 2 - lu'q 2 - p' 2 ) , (64) 

where Hb denotes the bulk Hamiltonian, namely the non-derivative terms in Eq. ()62j) . Thus we 
can employ Eq. Q41|l to calculate the surface tension here also. However, it is not possible to reduce 
the calculation of the surface energy to a simple quadrature over the density since here derivatives 
of three fields are involved. Note that in the relativistic case the Q's are the coefficients of the field 
gradients with 

Qaa = ~Qu>u) = ~Qpp = 1 j (65) 

all other coefficients being zero. 



B. Relationship to the Liquid Droplet Model of Nuclei 



The empirical liquid droplet approach, originally formulated in Ref. [22j|, provides a useful 
framework for the description of nuclei in the equation of state relevant for astrophysical simulations 
of supernovae and neutron stars. The droplet approach for isolated nuclei can be easily extended to 
the case in which nuclei are immersed in a dense medium comprised of electrons, positrons, protons, 
neutrons and alpha particles [3]. Historically, the parameters of the droplet model have been 
established from fits to binding energy data for laboratory nuclei (cf., Ref. |29|). Although such 
fits cannot independently establish reliable values for the volume and surface symmetry coefficients, 
they can determine a strong correlation between them. This correlation, however, depends on the 
formulation of the droplet model. 

The simplest droplet model consists of an extension of the Bethe-von Weizsacker liquid drop 
model to incorporate the surface asymmetry. As demonstrated by Myers and Swiatecki j^, the 
nuclear energy can be written as 

z 2 z 4 / 3 

+E dif — + E ex -^ + aAA^ 2 . (66) 

In this expression, B ~ 16 MeV is the binding energy per particle of bulk symmetric matter 
at saturation, E s ,Ec,Edif and E ex are coefficients for the surface energy of symmetric matter, 
the Coulomb energy of a uniformly charged sphere, the diffuseness correction and the exchange 
correction to the Coulomb energy, respectively. The last term represents pairing corrections, where 
A is a constant and a = +1 for odd-odd nuclei, for odd-even nuclei, and —1 for even-even nuclei. 
The effects of curvature and higher-order terms are neglected. The quantity S* s is related to S s as 
described below. 

Care must be taken in making the assumption that the surface energy scales as A 2 / 3 and the 
Coulomb energy as Z 2 A~ 1 / 3 . Physically, the Coulomb energy scales as Z 2 /R p . The surface energy 
scales with the surface area, but its definition depends on the exact radius of the surface and is 
therefore ambiguous [28j |. Nevertheless, the total of the bulk plus surface energies is unambiguous. 
This ambiguous splitting of the bulk and surface energies can be resolved if the surface thermo- 
dynamic potential is used instead of the surface energy. In what follows, we describe two ways in 
which this can be done. 
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"fi n " Approach 



In this approach developed by Lattimer et al. |2g], R p is used as the reference surface. This 
choice was made since the Coulomb energy is naturally expressed in terms of this radius. Since, 
in general, R n > R p , the bulk energy of an additional N s neutrons in the neutron skin has to be 
included. The diffuseness and Coulomb exchange corrections are also naturally expressed in terms 
of R p . These definitions allow the mass formula to be simply written in terms of the variables A 
and 1=1 — 2Z/A. However, the connection between R p and A depends on the density (assumed 
to be uniform and equal to til) of the bulk interior fluid. This results in the modified surface 
symmetry energy parameter 

where the true saturation density uq is used in place of the fiducial density of 0.16 fm~ 3 . Different 
parameterizations of the nuclear force predict different values for the saturation density, and con- 
sequently for the interior densities Ul in nuclei. Using R to refer to R p in this approach, the total 
droplet energy becomes 



E(A,Z) = {-B + S v 5 2 L )(A-N s )+A7rR 2 a(pn)+^nN : 
3Z 2 e 2 7r 2 Z 2 e 2 d 2 3Z 4 / 3 e 2 / 3 \ 



> \ ^ . - 1 /9 / \ 

+ W? +aAA ' 112 - < 68 » 

where explicit expressions for Eq, E^if and E ex have been substituted. The quantity d ~ 0.55 fm is 
the surface diffuseness parameter. Note that the bulk energy per baryon in the interior, expanded 
to second order in 5l, is —B + S V 5 2 L , and that the total surface energy is 4irR 2 a(iif l ) + f^nN s , where 
the surface thermodynamic potential per unit area a is a function of the surface neutron chemical 
potential fj,^. Technically, the surface tension is a function of both \x n and \x P) but if Coulomb and 
other finite-size effects are neglected, the matter pressure vanishes throughout the surface and (i n 
and n p can each be expressed as a unique function of 5l. Thus, the surface tension in this limit can 
be expressed as a function of only one chemical potential. Since by definition only neutrons exist in 
the neutron skin, it is appropriate to take the surface tension as a function of \x n . Minimizing the 
total energy with respect to [i s n requires that fi n be the same in the surface as in the bulk interior, 
and its value is thus given by Eq. (|43|) . The function <r(^^) is determined through the auxiliary 
variable 5l using a = do — The number of surface neutrons follows from minimizing the total 

energy with respect to N s and leads to the thermodynamic identity 

N s = -AnR 2 ^ = 4^ 2 -^— . (69) 

Noting that the local asymmetry in the bulk interior is 

A-N 8 -2Z 

6l = A-N s ' (70) 

one has that N s = A(I — #l)/(1 — 5l), where the global asymmetry is / = (N — Z)/A. Thus, it is 
possible to establish that 



S L =I 



4vr R 2 a s 
l H - 

S V A 



(71) 
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In order to express the modified liquid droplet energy in terms of A and /, it is necessary to 
relate the nuclear proton radius R to A, the interior density til and the asymmetry /. This is 
accomplished with the implicit relation 



—R 3 n L = A-N S 



A ^k= A{1 - !) 



i 



AI 



A + 4ttR 2 (ts/S v 



(72) 



This is a quintic equation for R. It is convenient to write its solution in terms of a function 
v(A,I,til) such that 



R = r [Av(A,I,n L )(0.16 fin" 



n 



-lil/3 



(73) 



where, as before, (47T7q/3)(0.16 fm~ 3 ) = 1. We also define S s = ^r^a^ to provide a measure of 
the surface symmetry parameter that is independent of the interior density, unlike S* in the simple 
liquid droplet model. In the case of symmetric matter, the quantity v(A, I, til) would simply equal 
unity; for asymmetric matter, it also contains the effects of the neutron skin. The implicit equation 
for v becomes 



(v - 1)(1 + av 2/3 — I) = -alv 2/3 



where 



a 



S„AV* 



0.16 fin" 



2/3 



(74) 



(75) 



In practice, solving the quintic equation for v can be avoided by employing its expansion 



1 



a 



a 



l + a 



(l + a) ; 



1 + ia] 



a 



1 2 2 2 1 
1 + -a + -a 2 + -a* P + 
3 3 9 



(76) 



in terms of the small quantity /. Eliminating N s ,5l and R, and setting w = t>(0.16 fm 3 )/n^, the 
liquid droplet mass formula becomes 



E(A, Z)= - BA + 4TTrla (wA) 2/3 + S V AI 2 



1 + 



S<w 2 / 3 



i -i 



S V A^ 



+ 



3e 2 Z 2 



(-) 



1/3 



2tqwA 



4r 



2-k J \wAJ 



aAA~ 1 / 2 . 



(77) 



For a given A and Z, therefore, the modified liquid droplet energy is a function of the six parameters 
B, o"o, S v , S s , til and A. The neutron skin thickness is given by 



2 r w S s 6 L 

3 v S v l-5 2 L 



with 5l = I 



1 + 



S s w 2 / 3 
S V AV3 



i -i 



(78) 



A shortcoming of the droplet model described above is it employs the surface tension along the 
coexistence curve for which the pressure is zero and Coulomb effects are ignored. Thus, in a nucleus 
with N = Z, it predicts that there should be no neutron skin. However, the Coulomb repulsion 
of protons results in a lowering of the proton density inside the nucleus and the develo pm ent 
of a proton skin. The following section describes the approach adopted by Danielewicz |33j to 
incorporate this Coulomb effect in the droplet model. 
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2. "/j, a " Approach 

In this approach, the surface tension is parameterized in terms of [i a = (fi n — n p )/2 instead of 
\x n . The quantity [i a is antisymmetric in neutrons and protons. In this way, the physical feature 
that the proton chemical potentials in the two bulk phases across a boundary differ is captured. 
Since neutrons and protons are treated symmetrically, the surface term n^Ng in the mass formula 
(|68|) is replaced by 

(/4 + fi)(N. + Z s )/2 + Ox- - ^)(JV a - Z.)/2 . (79) 

The surface thermodynamic potential is to be treated as function of fi s a alone. A result of these 
definitions is that minimization of the total energy with respect to fi s n — /i p and \i s n + [i s results in 
the surface asymmetry density being proportional to N s — Z s oc da/d^, s a , whereas the total surface 
density is proportional to N s + Z s oc da/dfi s and will vanish. Thus, Z s = —N s in this approach. 
In addition, the nuclear radius now satisfies 47rR s ni/3 = A. 

The possibility that Z s protons exist in the surface must be taken into account in the Coulomb 
energy. Then the Coulomb energy becomes 



to which diffuseness and Coulomb exchange corrections must be added. We note that Danielewicz 
kept an additional term in the expansion in Z s /Z, but this has only a marginal effect on the 
results. Minimizing the total energy 



E(A,Z)= ( -B + SJi){A -N s -Z s )+4irR 2 a(n s a ) + fi s (N s + Z s ) + fi s a (N s -Z s 
3Z 2 e 2 ( Z s \ 7T 2 Z 2 e 2 d 2 3Z 4 / 3 e 2 / 3 \ 2 / 3 A A 1/2 



(81) 



with respect to N s — Z s results in 



Using 



/4 = 2S^-^§. (82) 



a = a - a s (n s a /2S v ) 2 , (83) 

one finds that 

N s -Z s = 47rR 2 as^ = aA(5 L - (3) 

= A(I - 5 L ) = ocA 1 -^ , (84) 
I + aP 

fa = -j— — , (85) 
1 + a 

where a = / L7rR 2 as/(S v A) has the same definition as in Eq. (|75|) and f3 = e 2 Z/(20RS v ). The 
neutron skin thickness is 



2 N s -Z s = 2r S s 5 L -p 

3 A(l-S 2 ) 3u S v l-S 2 



t = R n -R p = -R-^ = -— Tg-, (86) 
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where uq = n^/0.16 fm 3 . This shows explicitly how the Coulomb repulsion reduces the skin 
thickness. Finally, the total energy can be written as 



E(A.Z) = -i M + 4, r M-f 3 + ^ /2 + a/J(2/ - /5) 



+ 



3e 2 Z 2 /w \ 1/3 



5r V A- 



UqJ 1 + a 

5n 2 d 2 /uo\ 2 / 3 5 ( 3 ^ 2 / 3 



6r 2 \AJ 4\2irZj 



+ aA.4" 1 / 2 . (87) 



IV. RESULTS 
A. Selection of Models and Their Parameters 

In this section, we wish to establish some generic trends that emerge from calculations based 
on the potential model Hamiltonian in Eq. © and the field theoretical Lagrangian in Eq. (|13j) . 
Toward this end, the various coupling strengths that enter in these two approaches are chosen so 
as to reproduce the empirical properties of 

equilibrium binding energy : B = — 16 ± 1 MeV , 

equilibrium density : uq = 0.16 ± 0.01 fm -3 , 

incompressibility : K = (200 - 300) MeV , 

Landau effective mass : m* L = (0.6 — 1.0) M , and 

symmetry energy : S v = (25 — 35) MeV . (88) 

In the case of potential models, calculations of nuclei are performed using the Hartree-Fock- 
Bogoliubov approach 21 ] that includes pairing interactions. Hartree calculations are employed 



in the field theoretical approach, as a treatment of the exchange (Fock) terms is considerably more 
complicated than in the potential model approach. In both approaches, we require that the binding 
energy and the charge radii of closed-shell nuclei are reproduced to within 2% of the measured 
values. The scalar meson mass in the field theoretical approach was restricted to lie between 450 
and 550 MeV. Insofar as fits to the valence single particle energies of closed shell nuclei are known 
to require a detailed treatment of correlations beyond the mean field level (such as short range and 
RPA correlations that are not considered here) we have chosen to tolerate slight deviations 
from the measured energies. We note, however, that Todd and Piekarewicz |3(| have discussed 
a correlation between the neutron radius of lead and the binding energy of the valence orbitals. 
They found that smaller neutron radii lead to last occupied neutron orbitals that are more weakly 
bound, and therefore smaller neutron drip densities. 

In all cases considered the supranuclear EOS was constrained to yield a maximum neutron 
star mass of at least 1.44Mq, the larger of the accurately measured neutron star masses in the 



double neutron star binary PSR1913+16 (see Ref. |37( for a compilation of known masses). This 
constraint resulted in the elimination of several non-relativistic models in which either the pressure 
decreased as a function of increasing density (e.g. SkP jH^j) or the proton (and hence the electron) 
fraction vanished at a finite supranuclear density (e.g. SkT3 H^). In the latter case, the condition 
of beta-equilibrium requires positrons and anti-protons which potential models cannot account for 
in a natural manner. Recently, a systematic study of potential models has been performed by 
Stone et al. who adopt a similar procedure to select among 87 Skyrme parameterizations. In 
addition to imposing the constraints in Ref. [i^l for potential models, we have also removed all 
models that do not have reasonable values for the Landau parameters at saturation density j^jj. 
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Since few relativistic models meet our criteria, we have generated a new class of relativistic 
models, es25, es275, es30, es325 and es35 with the Lagrangian of Eq. (j!3j) . In these models only 
a2 and b\ in Eq. (|15|) were allowed to be non-zero. They were designed to meet our restrictions 
while offering some variation in the value of the symmetry energy. We have also exploited the 
full freedom of the function / in Eq. (|15|) to produce models SRI, SR2 and SR3 for which the 
symmetry energy has a substantially weaker density dependence than in typical relativistic models. 
This leads to a smaller skin thickness. Note that es25 and es275 also have smaller skin thicknesses 
than typical. 

A list of all of the models used in this work together with their saturation properties, surface 
energies, and skin thicknesses is given in Table in Appendix A. In Appendix B the coupling 
strengths of models SR2, es25, es30, and es35 are listed in Tables IVl | IVIII fVIIII and lLX| respectively. 

Results for the EOS of APR 

We present here results for nuclei and semi-infinite matter obtained with the potential (NRAPR) 
and field-theoretical (RAPR) model fits to the EOS of APR. These models satisfy all of the con- 
straints mentioned above. The values obtained for the binding energy /particle and the charge radii 
of closed-shell nuclei are given in Table IIVI The results of the two models differ by at most a few 
percent and compare quite well with the data. This agreement is gratifying insofar as the EOS of 
APR was obtained from many-body theory without reference to nuclei. The experimental situation 
regarding skin thickness is quite model dependent and we show a few representative values in Table 
IIVI In agreement with the prediction, most of the recent extracted values for the skin thickness in 
40 Ca are small and negative (however positive values were obtained in earlier work, see the review 
of Batty et. al. For 90 Zr the skin thickness is predicted to be 0.084, on average, which is 

close to the central experimental value. In 208 Pb the theoretical value ~ 0.2 is at the upper end of 
current experimental estimates. 



TABLE IV: Comparison of results from the potential and field theoretical approximations to the APR 
equation of state with experimental data for selected closed-shell nuclei. 



Nucleus 


Property 


Experiment 


Potential 


Field-theoretical 


208p b 


Charge radius (fm) 


5.50 [43] 


5.41 


5.41 




Binding energy per particle (McV) 


7.87 [44] 


7.87 


7.77 




Skin thickness (fm) 


0.12 ± 0.05, 0.20 ±0.04 
[10,45] 


0.19 


0.20 


90 Zr 


Charge radius (fm) 


4.27 [43] 


4.18 


4.17 




Binding energy per particle (McV) 


8.71 [44J 


8.88 


8.65 




Skin thickness (fm) 


0.09 ±0.07 [46] 


0.075 


0.093 


40 Ca 


Charge radius (fm) 


3.48 [43J 


3.40 


3.34 




Binding energy per particle (MeV) 


8.45 [44] 


8.89 


8.61 




Skin thickness (fm) 


-0.06 ± 0.05, -0.05 ±0.04 

[10,42 


-0.044 


-0.046 



In Fig. 0]the neutron and proton density distributions for 208 Pb calculated with the NRAPR 
and RAPR models are seen to agree reasonably well with each other. These results are compared 
with the corresponding results of semi-infinite matter calculations which also agree reasonably 
well amongst themselves. There is, however, a noticeable enhancement of the relativistic densities 
compared to the non-relativistic ones at the onset of the surface region (r ~ 5.3 fm). In order to 
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FIG. 4: Calculated neutron and proton density distributions for 208 Pb and for semi-infinite matter, 
inset compares the calculated charge distributions with data ■ 



The 



specify Sl for the semi- infinite matter calculation, the neutron and proton densities were averaged 
from the center of the lead nucleus (for each case) out to 2/3 of its rms charge radius. This procedure 
smooths the oscillations in the calculated den sity distributions and roughly takes into account the 
effects of short range and RPA correlations [3J|. Note that semi-infinite matter calculations in 
the Thomas Fermi approximation are inherently unreliable in the tail of the distribution. The 
inset shows that the comparison between the calculated charge distributions and the data is 
comparable to that of other calculations performed at the Hartree-Fock-Bogoliubov (for potential 
models) or Hartree (for field-theoretical models) mean field level 0,0. 

Figure [5] shows neutron and proton density profiles from semi-infinite matter calculations for 
a range of <5l's for the NRAPR and RAPR fits to the EOS of APR. The isospin dependence of 
the profiles generally matches well with that of the semi-infinite matter profiles calculated directly 
from the EOS of APR. Values of 5l larger than about 0.3 (not considered for the RAPR case) 
necessitate the presence of dripped neutrons since the neutron chemical potential is positive. Hence 
the neutron density remains finite as z 



oo. 
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FIG. 5: The semi-infinite matter neutron (solid curves) and proton (dashed curves) density profiles for the 
indicated values of the asymptotic isospin asymmetry, Sl, as a function of distance z. The profiles are 
invariant upon translations along the axis and have been suitably shifted for ease of display. From left to 
right, results are for the EOS's of the APR, NRAPR, and RAPR models. 




FIG. 6: Left: The ratio of surface to volume contributions, S s /S v , to the total symmetry energy calculated 
from finite nuclei versus that calculated from semi-infinite matter. Right: The ratio S s /S v calculated from 
the semi-analytic approximation, Eq. I|52[l. versus that obtained from semi-infinite matter using Eq. 1)41(1 . 
In both panels, the downward (upward) pointing triangle is the result for the NRAPR (RAPR) model. 
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FIG. 7: Correlation between S s /S v and SR derived from nuclei (left panel) and between S s /S v and y3/5 £ 
derived from semi-infinite nuclear matter (right panel). 



B. Semi-infinite matter and finite nuclei 

Figure H3 compares different methods of calculating the ratio of the surface to the volume sym- 
metry energy, S s /S v . The left panel shows the values calculated from nuclei using Eq. (|78|) (with 
5l specified by the averaging procedure described earlier) plotted against semi-infinite matter re- 
sults from Eqs. (|^T|) and (|35|) . The latter values are expected to be the most reliable. The right 
panel shows the semi-infinite matter results from the approximate relation in Eq. ()52j) (using, 
however, exact expressions for the Hamiltonian and the symmetry energy) plotted against the 
exact semi-infinite results described above. For field-theoretical models, the factor Qnn ~\~ Qnp 
needed in Eq. ()52|) is calculated from Eq. (|46|) . In this figure, and in the following figures, potential 
(field-theoretical) model results are denoted by filled squares (circles). The circled triangles give 
the APR results from the two fits. There is some scatter about the straight line that indicates 
perfect agreement between these methods, but the correlation is sufficiently good that it validates 
the expressions we have used for this ratio. Notice that, apart from the newly constructed SRI, 
SR2, and SR3 models, the field-theoretical models tend to have a larger value of S s /S v than the 
non-relativistic models. 

In Fig. [7| we display the correlation between S s /S v and the neutron skin thickness 5R, or 
equivalently \/3/5 t. In the left panel these quantities are derived from nuclei, whereas for the 
right panel they are obtained from semi-infinite matter calculations. For the latter, the definition 
of t in Eq. (jSHJ is employed with 5l obtained from the lead nucleus (see above). Both panels 
display a good linear correlation, as suggested by Eq. and the aforementioned tendency of 

field-theoretical models to have larger values of S s /S v translates into a tendency toward larger 
values for the skin thickness. 
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C. Correlations and Their Origins 

In this section, we analyze in some detail the correlations that exist between quantities that are 
accessible through either laboratory experiments or astronomical observations that shed light on 
isospin-dependent interactions in nucleonic matter. Our emphasis will be on uncovering the causes 
of these correlations. 



1. Nuclear masses and the correlation between S s and S, 




20 25 30 35 20 25 30 35 

S v (MeV) S v (MeV) 



FIG. 8: Correlations in the symmetry energy parameters from liquid droplet fits to nuclear binding energies. 
The left (right) panel is for the '/i Q ' ('/in') method. Filled-in contours refer to the droplet models that have 
average errors 0.01, 0.02 and 0.03 McV/baryon, respectively, above the minimum values (marked by plus 
signs). The thick solid (dashed) line lies on the minimum x 2 valley for our (Daniclcwicz's) droplet models. 
Contours of skin thicknesses (SR) are shown for our droplet models; corresponding values arc indicated on 
Danielewicz's correlation with solid dots. The open circles indicate the minimum point found by Mpller et 
al. 0. 

The formalism for the modified liquid droplet models was discussed in Subsec. Ill B. Each 
model has six parameters and we wish to study the correlation between the volume and surface 
symmetry parameters, S v and S s (see also Lattimer and Swesty jHOl). Therefore as a function of 
S v and S s , we have minimized the quantity 

n 

X a = E(^-^) 2 /n, (89) 
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with respect to the other four parameters. Here Ef and Ei are the data |4J] and model energies 
of the ith nucleus and n = 2841 is the number of nuclei included in the fit (only experimentally 
determined masses with A > 20 are included; the results for a smaller set of nuclei with A > 40 
are not significantly different). In either case, a nearly linear behavior is found for the contour of 
minimum x 2 between S s /S v and S v , as shown in Fig. |HJ The equations describing the correlations 
are given by 

S s /S v = -5.253 + 0.2545,, for fj, a method , 

S s /S v = -3.453 + 0.1635^ for /j, n method . (90) 

The best-fit values of S v and S s /S v are about 27.3 (24.1) MeV and 1.68 (0.545), respectively, for 
the £ /i a ' ('/i n ') methods for the modified droplet models. 

A relevant question is the range over which values of S v and S s provide fits that differ by a 
statistically insignificant amount from the best fit. Traditional mass formula fits were developed to 
assist in interpolating within or to extrapolate to nearby regimes outside the ranges of measured 
nuclear masses. The experimental errors of measured masses are of order keV's per nucleus, so 
interpolations to the same accuracy are desirable. However, for many astrophysical applications, 
which involve extrapolations to ranges of A and Z well beyond the measured ranges, a model 
accuracy of a few hundredths of an MeV/baryon would be acceptable. Such an error in the 
computed energy per particle of a nucleus would translate into a similar error in the neutron or 
proton chemical potentials. Errors in chemical potentials generated by the use of the "single- 
nucleus" approximation, an essential assumption of liquid droplet high-density EOS models, are of 
order T/A [5l|, or several hundredths of an MeV, for example. 

At the best fit for the droplet models, for which the value of x 2 is Xo — 6 ; the mean difference, 
or error, between the experimental and predicted nuclear energy per baryon is 5 ~ 0.028 MeV, 
where 5 is defined by 

n 

xMMA) 2 /n, (91) 

i 

and ri is a random number between —1 and 1. Therefore, increasing the mean error by a further 
amount e per baryon would increase x 2 to a value 



.2 



X 



J2[r l (5 + e)A i ] 2 /n. (92) 



Contours of mean errors e of 0.01, 0.02 and 0.03 MeV per baryon above the minimum are shown in 
Fig. |H1 For the '/%' modified droplet model, and for an additional mean error of 0.01 MeV/baryon, 
the allowable excursions in S v are about (+5,-4) MeV, and in S s /S v , about (+0.9,-0.5), even 
though the excursions must be constrained by the obvious correlation. However, for the model, 
the allowed excursions in parameters are much larger. It is apparent that fitting nuclear energies 
alone cannot constrain these parameters, even though a relatively tight correlation between them 
exists. 

The origin of the correlations in Fig. |H1 can be most easily understood by examining the func- 
tional form of the total symmetry energy in the simple droplet model: 

-l 



E S ym(A, I) — S V AI 



1 + 



(93) 



Assuming that a linear correlation Sg/S v = a + bS v exists between the parameters, minimization 
of Eq. (|93|) with respect to S v yields 

A-^Y 1 = -5.31 , (94) 
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where the average is mass-weighted over the entire range of nuclear masses used in this study. 
Noting that 

= I ( Mn _ „ p + 2EcZA~^ + 2E dlf ZA~ 1 + \E ex Z^ A~^) , (95) 
where « n — fi p can be taken directly from nuclear masses as 

fin - Mp - + 1, Z) - E(N - 1, Z) - E(Z + 1, iV) + £(Z - 1, iV)] , (96) 

one can estimate the average symmetry energy, from which it follows that 

b = (/ 2 ) (E^y 1 (^A-^y 1 ~ 0.223 . (97) 

Thus, one has almost precisely the linear relation shown in Eq. (|9()[) for th e 'u „' droplet correlation. 

A method similar to the '/Xq,' method was suggested by Danielewicz [38|. His results resemble 
ours, but slight differences in the formulation exist, the most notable being that we account for 
variation of the central nuclear density in the fitting (i.e., the ratio n^/0.16 fm~ 3 = no in Eqs. 
m , m and m). * addition, Ret fitted parameters with an absolute energy deviation 
minimization, rather than the minimization of x 2 that we employ. Danielewicz determined best-fit 
values of S s /S v ~ 1.7 and S v ~ 27 MeV and found a correlation slope of 0.32 (compared to 0.25 
in our study); Danielewicz's fit is indicated by the dashed lines in Fig. |HJ We have reproduced 
Danielewicz's model fits by employing his stated procedures. It is evident that the position of 
the minimum and the slope of the correlation are affected by differences in the fitting procedure 
and the droplet formulae. Nevertheless, to within an average error of 0.005 MeV/baryon, the 
differences found are not significant. Importantly, Ref. found, as do we, that large excursions 
in parameter space are not excluded within the droplet model. The range of allowed excursions 
are much smaller for the modified droplet model. 

It is interesting that the model of Moller et al. gives S v = 30.8 MeV and S s /S v = 1.62. 
This point is indicated by a small open circle in the left and right panels of Fig. |H1 Clearly it 
lies much closer to the droplet correlation than the '/i a ' correlation. Although this point is 
displaced from the position of the best fit for the '/x n ' correlation by over 6 MeV in S v , its average 
error differs from the best fit by less than 0.02 MeV/baryon. In contrast, this point differs by about 
0.04 MeV/baryon from the best fit established with the '/j a ' correlation. We note that the model 
of Ref. |49J contains several additional effects compared to our modified droplet model, including 
terms representing compression, curvature, deformation and shell effects. It is thus not surprising 
that the best-fit parameters differ. However, in spite of the fact that the '/x a ' approach contains 
Coulomb corrections to the surface tension while the approach does not, the notable feature 
is that the model of M0ller et al. seems more consistent with the model than the l fi a ' model. 

We also compare the volume and surface symmetry coefficients determined for the potential and 
field-theoretical models developed in this paper with those of the published models utilized here and 
with the results of mass formula fits to nuclear energies. The parameters of the models developed 
in this paper have been fit to only four closed shell nuclei, while many of the previously-published 
forces have been fit to properties of a wider range of nuclei that span a larger range of asymmetries 
and masses. Thus, we anticipate that a correlation established from mass fits should be apparent 
when comparing previously-published forces, but not necessarily by the forces developed in this 
paper. This is precisely what is found, as indicated in Fig. |§J It is important to note that most 
of the previously-published models follow the correlation of the modified droplet model as 
opposed to the steeper '/i Q ' droplet correlation. This is in spite of the fact that the '/j a ' model 
takes into account the polarization of the nuclear interior due to Coulomb forces which the 
model ignores. 
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FIG. 9: The ratio of the surface and volume coefficients of the symmetry energy versus the volume coefficient 
for the models used in this work. The closed symbols show the exact values obtained from a quadratic fit 
to the surface tension in semi-infinite nuclear matter, whereas the open symbols use the value of S s /S v 
obtained via Eq. (|52fl . The solid (dot-dashed) line is the linear correlation from the '/i„' Cm™') droplet 
model Eq. (|90|) and the short-dashed line is the linear correlation from the droplet model of Ref. ■ 



These results highlight a major difference between the traditional use of droplet mass formula 
fits and that of current astrophysical applications. Historically, droplet mass formulae were used to 
interpolate within known nuclei to establish relatively precise values. In astrophysics, however, it is 
necessary to extrapolate to extreme conditions of asymmetry. Extrapolations depend significantly 
upon the correlation between S s /S v and S v . Given that the EOS's in current applications to 
supernova and neutron star matter are based on a modified droplet approach, it is necessary to 
use consistent parameter sets obeying the proper correlation. Otherwise, connections between 
quantities such as the neutron skin thickness and the neutron star radius could be incorrectly 
interpreted. 

Fig. El also displays the neutron skin thickness SR as a function of S s and S v for the two 
approaches. In the '// n ' approach, if the dependence upon no were to be ignored, contours of fixed 
5R would appear as horizontal lines since 5R would be a function of S s /S v alone. This would also 
be the case in the approach, if the dependence upon both uo and the Coulomb term (3 were 
ignored. The dependence of 5R on (3 is such that d(5R)/dS v < if the uo dependence is ignored, 
implying a positive slope for the 5R contours. Indeed, we find the slope of the 5R contours is larger 
in the approach than in the '/i n ' approach. However, Ref. finds that the 5R contours have 
negative slope; this discrepancy remains to be resolved. 

The relationship between the skin thickness and the ratio S s /S v was studied numerically using 
a variational approach by Bodmer jH^]. The results of his work can be understood from our 
analytical analysis in Section III C. Ref. [H^ also emphasizes the interplay between the symmetry 
and Coulomb energies. Ignoring the latter, as in the '/i n ' method, leads to larger values of 5R. 
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FIG. 10: The neutron skin thickness SR of finite nuclei versus the pressure of /^-equilibrated matter at a 
density of 0.1 fm -3 . 



2. The neutron skin thickness versus the pressure of subnuclear neutron-star matter 

Typel and Brown 0, Q| have noticed a correlation between the skin thickness and the pressure 
of pure neutron matter at a density of n = 0.1 fm -3 . To the extent that this correlation can be 
applied, a measurement of 5R will help to establish an empirical calibration point for the pressure 
of neutron star matter at subnuclear densities. Coupled with a neutron star radius measurement 
(see the discussion below), that could allow the inference of the pressure at supranuclear densities, 
the Typel-Brown correlation would be valuable in establishing the pressure-density relationship 
over a wide range of densities inside neutron stars. 

The connection between the neutron skin thickness and the symmetry energy has been known 
from Bodmer's work in the 60's pSl ]. The Typel-Brown correlation between 5R and P(5?io/8), 
which is closely related to the derivative of the symmetry energy, demonstrates clearly the new 
information that could be obtained by an accurate measurement of the skin thickness in 208 Pb. 
Typel and Brown considered a large set of Skyrme parameterizations including those that are 
inadequate to describe neutron-star matter |4(|. Such forces lead to small skin thicknesses. None 
of the field-theoretical models considered in Ref. 0] yields a skin thickness smaller than 0.2 fm. By 
using the function / in Eq. lfT3*|) . we have constructed examples of field-theoretical models, SRI, 
SR2, and SR3 which give skin thicknesses which are as small as those obtained from some potential 
models (~ 0.15 fm). Our models es25 and es275 also have small skin thicknesses. However, as 
pointed out by Furnstahl there are other possible ways in which this could be achieved, 
including the addition of a 5 meson or a more complete treatment of pions. If measurements in 
208 Pb give 5R < 0.15 fm, then Skyrme-like parameterizations have to be refined at high density to 
properly describe neutron stars. By the same token, field-theoretical models will need modifications 
to better describe low-density matter. 

We show in Fig. ^] the relation between the skin thickness 5R and the pressure of beta equili- 
brated matter at 0.1 fm -3 . We have verified that the correlation is similar also at a fixed fraction, 
5/8, of the saturation density, for a variety of potential and field-theoretical models. 
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Some insight into the origin of this correlation can be gained by examining how the beta- 
equilibrated pressure Pp (the abscissa in Fig. Ilfljl and the neutron skin 5R (the ordinate) vary with 
quantities that depend upon isospin. The pressure of cold, beta-stable neutron-star matter can be 
written as 

Pp(n,x) = n 2 [E'(n, 1/2) + E' sym (n)(l - 2x) 2 ] +P e + P M , (98) 

where x = n p /n is the proton fraction, E(n, 1/2) is the energy per particle of symmetric matter, 
E sy m(n) is the bulk symmetry energy (primes denote derivatives with respect to density), and the 
last two terms give leptonic (electron and muon) contributions. For beta-equilibrated matter, 

- Hp = -dE/dx ^ 4(1 - 2x)E sym {n) 

= He = Hfl » (") 

where /x e and are the chemical potentials of the electron and muon, respectively. This relation, 
together with the charge neutrality condition, n p = n e + n„, permits the evaluation of the equilib- 
rium proton fraction. Since muons begin to appear in matter for n > hq JHH], the only leptons at 
subnuclear densities are electrons. The electron pressure is 

P e = \n e He = \nxn e = nx(l - 2x)E sym (n) , (100) 

where we have used n e = n p for n < uq and Eq. I|99|l. Utilizing this relation in Eq. (|98|). we can 
write 



Pp{n, x) = P(n, 1/2) + n(l - 2x)E sym (n) 



{l _ 2x) d^E (n) +x 
dlnn 



(101) 



where the first term is the pressure of symmetric nuclear matter. Thus, the total pressure can 
be written at a particular density in terms of fundamental nuclear parameters. For densities at 
or below saturation density, the equilibrium proton fraction, x, ~ (AE sym /hc) 3 /(3ir 2 n), is very 
small so that the beta equilibrated pressure and the pure neutron matter pressure only differ by 
a small amount that depends upon the symmetry properties of the force. Thus, for n < uq, 
Pp{n) = P(n, 1/2) + n 2 E' sym , where the second term provides the bulk of the contribution. Thus 
an approximately linear correlation would result if in Fig. 5R were plotted against E' sym (n) at 
densities n = 0.1 fm~ 3 or n = 5no/8; such a plot is given in Ref. jo^ |. 

The dependence of 5R on the isospin asymmetry can be seen explicitly in Eq. (|78jl for nuclei and 
in Eq. (|54|) for the skin thickness of semi-infinite matter. The latter equation can be manipulated 
to read 

}u 1 / 2 [h B (u,5 = 0) + ZBf 2 [S v /E sym (u)-l]du 





where u = njn^ and Hb{u,5 = 0) = TCb(u,5 = 0)/n. Since t = y/5/3 5R, this relation illustrates 
that the neutron skin of nuclei is proportional to a specific average of [S v /E sym (u) — 1] in the 
nuclear surface, the averaging function involving the square root of the energy per baryon of 
isospin symmetric matter. 

Equations (|98|) and (|102|) show that the density dependence of the symmetry energy is the 
principal cause of the Typel-Brown correlation. We therefore display the symmetry energies for 
the models considered in this work in Fig. 1111 The thick line indicates the symmetry energy of 
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n (fm" 3 ) 

FIG. 11: The symmetry energy versus density for various equations of state. Solid (dashed) lines are for 
field-theoretical (potential) models. The thick solid line shows the APR symmetry energy for the low density 
phase. The square (circular) markers indicate the densities at which potential (ficld-thcorctical) models allow 
the direct Urea process to occur. All models are constrained as discussed in the text. The arrow identifies 
the NRAPR (dotted line) and RAPR (solid line) models which have nearly identical symmetry energies at 
n = 0.5 fm~ 3 . 

the APR equation of state for the low density phase (this is strictly only applicable for n < 0.25 
fm~ 3 for which region the NRAPR and RAPR fits show good agreement). Generally the symmetry 
energies of the field-theoretical models (solid lines) vary more rapidly with density than those of 
potential models (dashed lines). Three exceptions, models SRI, SR2, and SR3, were constructed 
in this work by a suitable choice of the function / in Eq. p5|). These field-theoretical models 
mimic the density dependence of potential models near the saturation density, but with somewhat 
different behavior at supranuclear densities. Notice that the requirements that we have placed on 
our models strongly constrain the symmetry energy at an average nuclear density n ~ 0.1 fm~ 3 , 
in agreement with Brown 0| and Horowitz and Piekarewicz Q (see also Ref. |i3|). 

For the models considered in Fig. 111! the variation of 5R with Pp shown in Fig. confirms the 
Typel-Brown correlation. Additional points that emerge from our analysis include: 

(1) At the densities of interest, the pressure difference between pure neutron matter and beta- 
equilibrated matter is small enough that the correlation is not significantly impacted. 

(2) In general, field-theoretical models yield 5R > 0.2 fm. The exceptions are models SRI, SR2, 
SR3, es25, and es275 that were designed to have special symmetry properties that lead to 5R < 0.2 
fm. 

(3) If the neutron skin thickness is measured with high accuracy, the pressure of neutron-star matter 
at this density could be determined to within (20-25)%. For a fiducial R n — R p = 0.2 ± 0.025 fm, 
this correlation predicts that Pp = (0.9 ± 0.3) MeV/fm 3 ; for a lower mean value of R n — R p , 
the predicted uncertainty in the pressure is larger. Employing a wider variety of non-relativistic 
models than those used by Typel and Brown results in an enlarged spread in results for models 
with relatively small values of 5R. We emphasize that the linear correlations observed in Fig. ^] 
are significantly worsened for slight excursions from the fiducial densities chosen. 

Furnstahl [541 ] has presented a correlation between the neutron form factor and the neutron 
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radius. He also demonstrated correlations between the skin thickness and, separately, the linear 
density dependence of the symmetry energy, the quadratic density dependence and also d 2 no/da 2 . 
These three quantities are all closely related to P(5no/8) as shown above. 

Implications for neutrino emission from neutron stars 

The density dependence of the symmetry energy also plays an important role in determining 
whether or not the simplest possible v emitting processes, the so-called direct Urea processes 

p + i — > n + U£ and n — > p + I + Vg, (103) 

where I is either an electron or a muon, occur in charge neutral neutron star matter. These direct 
Urea processes can occur whenever energy and momentum conservation is simultaneously satisfied 
among n, p and I (for temperatures of interest to long-term cooling, neutrinos would have left the 
star). In a mixture of neutrons, protons and electrons, the required proton fraction x for the direct 
Urea process to occur is 1/9. When this mixture is in beta equilibrium, x satisfies 

x ~ 0.048 (E sym (n)/S v f (n /n) (1 - 2xf , (104) 

which highlights the role of the density dependent symmetry energy. For example, for a linear 
symmetry energy, E sym {n) = S v (n/no), the onset of the direct Urea process occurs for n c /riQ ~ 2.2 
with S v = 30 MeV. In the presence of muons, which begin to appear around no, the required 
proton fraction is about 0.14 and hence the threshold density is somewhat larger. In the case that 
E sy m(n) increases less than linearly with density, the direct Urea threshold density is larger. If 
the Urea threshold density is less than the star's central density, the neutron star will rapidly cool 
because of large energy losses due to neutrino emission: the star's interior temperature T will drop 
below 10 9 K in minutes and reach 10 7 K in about a hundred years. This is the so-called rapid 
cooling paradigm. If the threshold density is not reached below the central density, or if the direct 
Urea process is suppressed due to nucleon superfluidity (calculated energy gaps are about an MeV 
or less), cooling instead proceeds through the significantly less rapid modified Urea processes 

n + (n,p) — > p + (n,p) + e~ + v e and p + (n,p) — > n + (n,p) + e + + u e , (105) 

in which an additional nucleon enables momentum conservation. As the temperature approaches 
the superfluid critical temperature, cooling also takes place through neutrino pair emission from 
the Cooper-pair breaking and formation processes which are more efficient than the modified 
Urea process, but significantly less efficient than the direct Urea process. If the density dependence 
of the nuclear symmetry energy is relatively weak, only the less rapid neutrino cooling processes 
can occur. 

The filled squares and circles in Fig. ^2 indicate the densities at which the direct Urea process is 
allowed by energy and momentum conservation conditions. Some of the models whose symmetry 
energy exhibits a moderate density dependence allow the direct Urea process to take place, but at 
densities beyond the maximum of 0.5 fm -3 plotted. For some other models the direct Urea process 
does not occur at all since it requires densities larger than the central densities of the maximum 
mass stars. 

Figure IT2*1 shows the skin thickness versus the threshold density for the direct Urea process (note 
that the NRAPR and RAPR skin thicknesses are plotted against n\j rca for the exact APR EOS.). 
The results in this figure supplement those of Horowitz and Piekarewicz by also including 
results of potential models. A determination that 5R has a value greater than 0.2 fm would imply 
that the nucleon direct Urea process should exist in 1.4 M neutron stars, unless it is suppressed 
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FIG. 12: Skin thickness versus the threshold density for the direct Urea process to occur in neutron stars. 

by superfluidity. For a recent assessment in this regard, see Ref. [Hi| which studied neutron star 
cooling in the absence of enhanced neutrino emission, and compared the results with data on 
cooling neutron stars. 



Implications for the crustal fraction of the moment of inertia 

Another astrophysical application concerns the thickness of the crust of a neutron star, which 
is related to the density dependence of the symmetry energy by physics similar to that of the 
thickness of the neutron skin of nuclei. The neutron star crust thickness might be measurable 
from observations of pulsar glitches, the occasional disruptions of the otherwise extremely regular 
pulsations from magnetized, rotating neutron stars. The canonical model of Link et al. [HO] suggests 
that glitches are due to the transfer of angular momentum from superfluid neutrons to normal 
matter in the neutron star crust, the region of the star containing nuclei and nucleons that have 
dripped out of nuclei. This region is bounded by the neutron drip density (~ 1.5 x 10~ 3 no — 4x 10 11 
g cm -3 ) and the transition density nt . f~ no/2 ~ 1.5 x 10 14 g cm" 3 ) at which nuclei merge into 
uniform nucleonic matter. Link et al. [60j concluded from glitches of the Vela pulsar that at least 
1.4% of the total moment of inertia resides in the crust of the Vela pulsar. The fraction of the 
star's moment of inertia contained in the crust is connected to the crust's thickness, and both are 
controlled by the pressure of matter, Pt, at the transition density as well as the neutron star's mass 
M and radius R. The fractional moment of inertia, AI/I, where Al is the moment of inertia in 
the star's crust and / is the star's total moment of inertia, can be expressed as j3| : 

Al _ 287rP t R 4 (1 - 1.67/3 - 0.6/3 2 ) 
T~ ~ 3GM 2 

where (3 = GM/Rc 2 . 

The core-crust transition density is not too far from 0.1 fm -3 , and matter at these densities is 
nearly pure neutron matter, which suggests that the Typel-Brown correlation between nt and 5R 
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FIG. 13: Neutron skin thickness versus the core-crust interface density nt (left panel) and the fractional 
moment of inertia SI /I in the crust of a 1.4Mq neutron star (right panel). The dashed line for AI/I = 0.014 
is deduced for the Vela pulsar, assuming a mass for Vela of 1.4M Q . 



has astrophysical significance. However, the Typel-Brown correlation refers to the fiducial density, 
0.1 fm" 3 . Thus, even if the transition density was always 0.1 fm 3 , the dependence of R (for a 
given M) on the symmetry energy has to be taken into account also. The major dependence on 
the EOS is therefore through the quantity P t R A /M 2 . 

A reasonable approximation to the core-crust boundary can be obtained by considering the zero 
temperature phase equilibrium between two homogenous nucleonic fluids, a dense phase approx- 
imating matter inside nuclei (labelled by i) and a dilute phase approximating the drip nucleons 
outside nuclei (labelled o). Finite-size effects, such as those due to the surface and Coulomb energy 
can be included as in Ref. [63], but are ignored here to obtain a first orientation. The electrons 
can be treated as a free Fermi gas which uniformly fills space in both phases and we define Y e to 
be the ratio of electron and baryon densities. 

Phase equilibrium follows from minimization of the free energy density allowing for the two 
conserved quantitites of baryon number and charge. We define u to be the fractional volume 
occupied by the dense phase i. The dense phase i has a proton fraction xf, the dilute phase o has 
a proton fraction x Q . Phase equilibrium, while maintaining beta equilibrium, is then described by 

n = riiu + n (l — u) , nY e = riiXiU + n x {l — u) , Pj = P Q , 

The transition density occurs at the baryon density n = nt for which u — > 0. 

Figure fTTH displays the relation between 5R and nt (left panel) and the relation between 5R and 
AI/I (right panel) for a 1.4 Mq star. The transition density varies over a range of about a factor of 
2 in density, from 0.05 fm~ 3 to 0.12 fm~ 3 , and there is a negative correlation between 5R and nt, but 
it is not very robust when all the models are considered. We note that Horowitz and Piekarewicz 0] 
computed the transition density as the density at which an instability to small perturbations occurs 
in a uniform liquid of nucleons and electrons (see also Ref. |25|). They obtained a much tighter 
correlation between nt and 5R by varying the symmetry energy in a few relativistic models. Using 
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FIG. 14: The quantity RP~ a as a function of the stellar radius R, for pressures P determined at 3/2, 2 
and 3 times equilibrium nuclear matter density. For each density, the best-fit value for the exponent a is as 
indicated. Results for 1.4Af Q (maximum mass) stars are in the left (right) panel. Circled stars indicate the 
results obtained with the APR equation of state. 

the exact APR EOS our method yields nt = 0.094 fm -3 in close agreement with Ref. |25| : this lies 
between the RAPR and NRAPR results plotted in Fig. ^3 being somewhat closer to the former. 

Although we do not display it, there is little correlation between P t and 5R. One would expect 
high values of S' v to correspond to high pressures in general. However, high values of S' v also tend to 
have low transition densities, so the effects somewhat cancel, leaving little net correlation between 
Pt and S^. However, as discussed above, the astrophysical observable AI/I is nearly proportional 
to the quantity PtR 4 /M 2 and not Pt alone. This is plotted in the right panel in Fig. 1131 Although 
the correlation is not very strong, it does suggest that, at least in the case of a 1.4M© Vela pulsar, 
a value of AI/I ~ 0.014 implies a skin thickness > 0.15 fm for 208 Pb. A larger mass for the Vela 
relaxes this lower bound toward lower values of 5R. 



3. The neutron star radius versus the pressure of supranuclear neutron-star matter 

Lattimer and Prakash found that the radius R of a neutron star exhibits a correlation that 
has the form of a power law: 

R ~ C(n, M) [P(n)]°- 23 -°- 26 , (108) 

where P(n) is the total pressure inclusive of leptonic contributions evaluated at a density n in the 
range 1 to 2nn, and C(n, M) is a number that depends on the density n at which the pressure was 
evaluated and on the stellar mass M. In Fig. El the left and right panels show this correlation as 
RP~ a versus R for stars of mass 1.4M© and M max , respectively, for the EOS's considered here and 
densities n = 1.5 — 3nn. In each case, the exponent a was determined by a least-squares analysis to 
give the best correlation. For the optimum exponent a = din R/ din P ~ 0.25, the radius increases 
very slowly with mass. As shown in Ref. 0], general relativity reduces the value of a from 1/2, the 
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FIG. 15: Left panel: The smallest radii of neutron stars given by our set of EOS's as a function of the 
corresponding maximum masses for field-theoretical and Skyrme models (respectively, solid and dashed 
lines labelled "max"). The radii for 1.4M Q neutron stars are labelled "1.4". The shaded regions indicate 
radii that are accessible for a particular class of EOS. Right panel: The mass-radius relation for the field- 
theoretical model SRI. The central densities of stars of different masses are as indicated. The inset shows 
the adiabatic index T versus density for this model. 



value characteristic of an n = 1 Newtonian polytrope. For the EOS of Buchdahl 62], this role of 
general relativity was demonstrated analytically 

In studying the Lattimer-Prakash radius-pressure correlation we have excluded cases in 
which the high density EOS was softened because of the presence of non-nucleonic degrees of 
freedom. The correlation, however, is not greatly improved by this omission, primarily because 
softening components do not often affect the equation of state strongly until the density exceeds 
two to three times the saturation density. 

The discussion in Sec. IIV C 21 showed that the pressure at densities close to equilibrium density 
is dominated by the derivative of the symmetry energy. Thus one could alternatively achieve a 
linear correlation by plotting E' sym (3no/2) versus i? 1//a . 

Sizes of neutron star radii in potential and field-theoretical models 

We have not considered in this work the possible presence of hyperons, Bose (pion or kaon) 
condensates or quarks. Such components generally soften the EOS which results in neutron stars 
with masses and radii that are smaller than those with nucleons-only matter. Our interest here lies 
in establishing the minimum radii that EOS's with only nucleonic matter yield. Fig. 1151 displays 
this minimum radius as a function of maximum mass for potential and field-theoretical EOS's that 
are constrained as described in Sec. IIV Al 

Each point on the solid line labeled "max" is the radius of the maximum mass configuration for 
the field-theoretical EOS which yields the smallest possible radius. The solid line labeled "1.4" is 
the radius of a 1.4Mq star for the same EOS. Analogous results for potential models are shown by 
the dashed lines. The hatched regions indicate where other reasonable potential or field-theoretical 
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FIG. 16: Calculated neutron skin thickness 8R of nuclei versus the radius of 1.4M© stars (left panel) and 
of maximum mass stars (right panel). The solid lines are described in the text. 

models may exist. For comparison, the APR equation of state supports a star of maximum mass 
2.2Mq for which the radius is 10.9 km. The field-theoretical models SRI, SR2, and SR3 have 
maximum masses of 1.44, 1.6, and 1.8 Mq, respectively. The corresponding radii are 

Ria = 9.84, 11.04, and 12.21 km for M = 1.4M 
Rmax = 9.17, 9.83, and 10.49 km for M = M max . (109) 

The mass-radius relation of the model SRI with a maximum mass of 1.44 Mq is shown in the right 
panel of Fig. 1151 Of the many field-theoretical models considered in this work, this model yields 
the smallest radius. 

4- The neutron skin thickness versus the neutron star radius 

We turn now to the Horowitz-Piekarewicz correlation for the models considered in this work. 
Our results are shown in Fig. ^] for a 1.4M Q (left panel) and for the maximum mass (right panel) 
configuration for each EOS. Note especially the results of models SRI, SR2, and SR3 which buck 
the trend of field-theoretical models that generally yield large neutron star radii. As noted in 
Ref. 0], large radii are a consequence of symmetry energies that vary rapidly with density beyond 
the nuclear saturation density. 

The Typel-Brown correlation demonstrates that the neutron skin thickness 5R is linearly cor- 
related with E' sym (5no/8), whereas the Lattimer-Prakash correlation for the neutron star radius 
establishes a linear correlation between R l l a and E' sym (3no/2). These two correlations can be 
synthesized by comparing E' sym at the two different densities. This is carried out in Fig. El 
where, notwithstanding some scatter in the results, there is an overall linear correlation between 
E' sym (5no / 8) and E' sym (3no/2). This suggests a relationship of the form 8R = a + bR}/ a . Deter- 
mining the coefficients by a least-squares fit to the results in Fig. Ilril yields 

SR = (0.09558 + 1.201 x 10" 5 R 1/0 - 281 ^) fm for M = 1.4M 
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FIG. 17: The derivative of the symmetry energy at 3no/2 versus that at 5no/8. 

5R = (0.1160 + 4.540 x 10~ 6 i? 1 / - 245 ) fm for M = M max , (110) 

These relationships are depicted by solid lines in Fig. Ilril For a fiducial 5R = 0.2 ± 0.025 fm, 
Ria = (13 ± 0.5) km and R m ax = (H ± 0.5) km. In the latter case, M max is found to vary up 
to 2.2M . These values give a measure of the extent to which the neutron star radius can be 
constrained if the neutron skin thickness is measured to the specified accuracy. For a lower mean 
value of 5R the uncertainty in the radius is larger. 

Schramm (6^1 further explored the Horowitz-Piekarewicz correlation 0,E3| between the neutron 
skin thickness and the neutron star radius by using covariant field-theoretical models based on 
SU(3) symmetry. He found that the correlation was weak because of the restricted variation 
allowed in the fourth-order vector-isovector couplings. In our work we have eliminated EOS's 
which do not satisfy minimal astrophysical and experimental nuclear constraints, which results in 
a tight correlation despite the enlarged freedom in the description of the symmetry energy. In fact, 
combining the correlations as discussed above, we find a skin thickness which scales approximately 
with the stellar radius raised to a power ~ 4. 



V. OTHER RELATED CORRELATIONS 

We turn now to address a few other correlations connected with isospin asymmetry. Supplements 
to our brief account here are contained in several reviews, some of which are alluded to below. 



A. Giant Dipole Resonances 

The analyses of giant resonances, particularly dipole resonances, have long served to delineate 
the role of volume and surface effects in nuclei; for a review see Ref. |65l |. In medium to heavy nuclei 
in which magnetic contributions are small, the inverse-energy-squared weighted photoabsorption 
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cross section can be related to the static polarizability p as 



a_ 2 = I du j = 2ir 2 (e 2 /hc)p . (Ill) 



In nuclei with mass number A > 100, cj_ 2 = (2.9 ± 0.2) A 5 / 3 fib MeV -1 parameterizes the data |6J, 
The dipole polarizability can be evaluated as the response of a nucleus to an external dipole 
field rfD, where if denotes the strength and D = (1/2) Ylt=i z i T t ls the dipole operator. Explicitly, 

\(0\D\n)\ 2 , . 

P = 2]T IV 1 1 n , (112) 

where \n) and u n are the eigenstates and eigenenergies of the nuclear Hamiltonian responsive to 
the dipole operator D. Microscopic RPA calculations [f3 of p employing Skyrme-like interactions 
reproduce the general trends of the data. The interplay between volume and surface effects are, 
however, difficult to extract from RPA calculations. Semiclassical methods, in which the energy 
density formalism is coupled with a hydrodynamic approach to describe collective excitations, 
have thus been employed to explore how the volume and surface symmetry energies, S v and S s , 
determine the value of p across the periodic table 0,0 0. Using the Hamiltonian density in 
Eq. (|T|). and writing the transition density as 5[p n (r) — p p (r)] = rf(f>(r) cos 9 , the polarizability is 
calculated from 

2-7T f 

P = y J Hr)r 3 dr , (113) 

where 4>(r) is the solution of the corresponding Euler-Lagrange (integro-differential) equation. It is 
found that similar results for p are obtained for correlated values of S v and S s . Numerical results 
for 40 Ca, 120 Sn, and 208 Pb show that the choice of S v in the range 27-42 MeV requires \S S /S V \ = 
1.2-2.2; lower values of S v correlate with lower values of S s for good fits The slope of this 

correlation can be compared to that found using ground state energetics in Sec. IV C 1 (see Fig. 
IHJ. It is in close agreement with the correlation obtained with the '/i n ' model. 

A qualitative understanding of this correlation can be gained by using a schematic symmetry 
energy density functional and a leptodermous expansion of the density [651 ] . whence one obtains 
the result 

-^Hl— ■)■ 

where (r 2 ) is the mean square radius of the nucleus and higher order terms contain corrections 
from the diffuseness and the skin thickness. Although values of S v = 32.5 MeV and \S S /S V \ = 2.2 
describe the data adequately [6f|, correlated variations in these numbers are allowed as pointed 
out in Ref. 

Sum rules have been particularly useful to relate experiments and theory in discussing the mean 



excitation energies, widths, and the spreading of the excitation strengths |65jl. The moments m„ 
of the strength function S(lu) = J2n>o \ ( n 1-^1 0)1 — w n ) defined by 

roc 

m p = S{uj)uj p duj = Y \(n\F\0)\ 2 ioP, (115) 

where F is the physical operator exciting the nucleus from its ground state |0) to its eigenstate \n), 
are especially helpful in this regard. For example, a good measure of the mean excitation energy is 
provided by E(D) = ^/m^JrnZx, for which results from RPA and hydrodynamic calculations have 



41 



been compared with data [65J]. Using a droplet model coupled with a hydrodynamic approach to 
excite the dipole resonance, Ref. [llJ] obtains 



E(D) 



6h 2 (l + K D ) S v 



where M is the nucleon mass. The quantity Kd is a model dependent enhancement factor char- 
acterizing the relative contribution of the nuclear interaction to the m\ sum rule and depends 
critically on the value of the energy up to which the energy integration is carried out in analyzing 
the experimental data. Using values of S v = 32.5 MeV and S s /S v ~ 2 from a droplet model 
fit to nuclear energies, and Kd = 0.2 corresponding to E max = 30 MeV, Ref. finds that a 
hydrodynamic approach is able to reproduce the experimental mean excitation energy of nuclei 
ranging from °Ca to 208 Pb reasonably well. Several Skyrme interactions with different values of 
S v and S s /S v are also able to account for the data to the same level of accuracy. However, since the 
values of S v and S s for these forces are correlated because they are fit to experimental masses, the 
additional constraint on the permitted ranges of S v and S s resulting from fitting dipole resonances 
does not seem to be significant. 

It must be emphasized that while a hydrodynamic approach is able to account for the gross 
features of nuclear mass dependence, several detailed features of the data such as strength frac- 
tionation, spreading widths, etc., are not naturally incorporated in its scope. For such details one 
must adopt a more microscopic approach that includes, for example, contributions from 2p-2h 
excitations, etc. For a detailed account, see, 

for example, Ref. H- 



B. Heavy-Ion Collisions 

1. Multi-fragmentation 

The breakup of excited nuclei into several smaller fragments during an intermediate-energy 
heavy-ion collision probes the phase diagram of nucleonic matter at sub-saturation density and 
moderate (~ 10 — 20 MeV) temperatures. In this region of the phase diagram the system is 
mechanically unstable if {dP / dn)T,x < 0, and/or chemically unstable if {d(i p / dx)x,p < 0. (A 
pedagogical account of such instabilities can be found in Ref. [z2]). These instabilities, which are 
directly related to the symmetry energy at sub-saturation densities j^l] , are believed to trigger the 
onset of multifragmentation. Because of the instabilities, matter separates into coexisting liquid 
and gas phases, which each have different proton fractions, i.e. "isospin fractionation" |7J]. This 
fractionation is observed in the isotopic yields which can potentially reveal information about the 
symmetry energy. Also, the scaling behavior of ratios of isotope yields measured in separate nuclear 
reactions, "isoscaling" , is sensitive to the symmetry energy |75L l76l| . This scaling is expressed in 
the empirically observed ratio of fragment yields from two similar systems with different neutron- 
to-proton ratios: 

Y 2 (N,Z)/Y 1 (N,Z) ~ exp aN+ ? z , (117) 

where the constants a and f3 can be related to the neutron and proton chemical potentials in 
a canonical ensemble description, and thus to the proton fraction of the source. To date, there 



are many suggestions of how the symmetry energy may affect multifragmentation |77(. Ongoing 
research is concerned with an extraction of reliable constraints on the symmetry energy from the 
presently available experimental information. 
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2. Collective Flow 

Nuclear collisions in the range Ei a \,/A = 0.5 — 2 GeV offer the possibility of pinning down the 
equation of state of matter above normal nuclear density (up to ~ 2 to 3no) from a study of matter, 
momentum, and energy flow of nucleons The observables confronted with theoretical analyses 
include (i) the mean transverse momentum per nucleon (p x )/A versus rapidity y/y pro j [zH- (h) 
flow angle from a sphericity analysis (iii) azimuthal distributions jsjj, and (iv) radial flow 
jj^ . Flow data gathered to date are largely for protons (as detection of neutrons is more difficult) 
and for collisions of laboratory nuclei in which the isospin asymmetry is not large. Theoretical 
calculations have generally been performed using Boltzmann-type kinetic equations. One such 
equation for the time evolution of the phase space distribution function f(f,p, t) of a nucleon that 
incorporates both the mean field U and a collision term with Pauli blocking of final states is (see, 
for example, Ref. 



^ + V P U ■ V r f - V r U -V p f= - J d 3 p 2 d 3 p 2 , cM^- v 12 (2tt) 3 5 3 (p + p 2 - p v - Pv) 

x [fh(i - h>)(i - h) - fvh>(i - /)(i - h)] ■ (us) 

Above, da NN /dn is the differential nucleon-nucleon cross-section and v 12 is the relative velocity. In 
general, the mean field U depends on both the density n and the momentum p of the nucleon. Equa- 
tion ()118j) contains effects due to both hard collisions and soft interactions, albeit at a semiclassical 
level. Theoretical studies that confronted data have thus far used isospin averaged nucleon-nucleon 
cross sections and mean fields of symmetric nuclear matter. It is now well established that much 
of the collective behavior observed in experiments stems from momentum dependent forces at play 
during the early stages of the collision jSJ]. The conclusion that has emerged from several studies 
is that as long as momentum dependent forces are employed in models that analyze the data, a 
symmetric matter compression modulus of ~ 220 MeV, as suggested by the analysis of the giant 
monopole resonance data [jjl ]. fits the heavy-ion data as well |8q . 

The prospects of rare isotope accelerators (RIA's) that can collide highly neutron-rich nuclei has 
spurred further work to study a system of neutrons and protons at high neutron excess 87 , 8^, 
Generalizing Eq. l|118j) to a mixture, the kinetic equation for neutrons is 



^ + V p U ■ V r f n - V r U ■ V p / n = J n = £ J ni , (119) 

i=n,p 

where J n describes collisions of a neutron with all other neutrons and protons. A similar equa- 
tion can be written down for protons with appropriate modifications. On the left hand side of 
each coupled equation the mean field U = U(n n ,n p ;p) depends explicitly on the neutron-proton 
asymmetry. The connection to the symmetry energy arises from the fact that U is obtained from 
a functional differentiation of a Hamiltonian density, such as that in Eq. Q. Examples of such 
mean fields may be found in Refs. [sli. Isol l9oT|. Observables that are expected to shed light on the 
influence of isospin asymmetry include neutron-proton differential flow and the ratio of free neutron 
to proton multiplicity as a function of transverse momentum at midrapidity. Experimental inves- 
tigations of these signatures await the development of RIA's at GeV energies. In this connection, 
it will be important to detect neutrons in addition to protons. 



VI. DISCUSSION AND CONCLUSIONS 



In this work, we have investigated the relationship between the symmetry energy of nucleonic 
matter and the neutron skin thicknesses of neutron-rich nuclei as well as the radii of neutron stars. 
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Precision measurements of neutron skins and neutron star radii are due to become available in 
the near future. We have studied how these measurements can constrain the dependence of the 
symmetry energy on baryon density, E sym (n), in the vicinity of the nuclear saturation density no- 
This knowledge is crucial to understanding many astrophysical phenomena, including neutron star 
evolution, supernova explosions and nucleosynthesis, and binary mergers involving neutron stars. 
Simulations of these phenomena involve extrapolations of the symmetry energy to supranuclear 
densities. 

A crucial aspect of our presentation is that care was taken to ensure that the nuclear force 
parameterizations, whether for potential or field-theoretical models, were constrained by fits to 
closed-shell nuclei. In addition, all models were constrained to yield a maximum neutron star mass 
of at least 1.44Mq, the larger of the accurately measured neutron star masses in the binary pulsar 
1913+16. Without such constraints, much weaker correlations relating the neutron skin thickness 
to astrophysical quantities such as neutron star radii would have resulted. 

Analytical representations of the relation between the symmetry energy and the neutron skin 
were developed for semi-infinite surfaces calculated in potential and field-theoretical models. To 
lowest order, the neutron skin thickness 5R = (r 2 ,) 1 ^ 2 — (r 2 ) 1 ^ 2 is proportional to 5lS s /S v where 
5l is the neutron excess in the nuclear center and S v [= E sym (no)] and S s are related to the 
expansion parameters of the volume and surface symmetry energies, respectively, in the liquid 
drop or droplet models of nuclei. This representation is validated by comparisons with results 
of finite nucleus calculations performed in the Hartree-Fock-Bogoliubov (for potential models) 
and Hartree (for field-theoretical models) approximations. For the first time, results for nuclei are 
presented for the Akmal-Pandharipande-Ravenhall EOS parameterized according to both potential 
and field-theoretical models. 

The semi- infinite surface representation predicts that the ratio S s /S v , or equivalently, the skin 
thickness 5R, is a particular average of the density-dependent symmetry energy in the surface. In 
addition, phenomenological comparisons indicate that there is a relatively tight correlation between 
5R and the pressure of neutron star matter at a typical average surface density of 5no/8. These 
results are not independent, but together imply that the determination of 5R offers a valuable 
constraint on the density dependence of the symmetry energy. It was explicitly shown that the 5R 
- pressure correlation is a consequence of a more general correlation between 5R and the derivative 
of the symmetry energy at the same density. In contrast to most, but not all, previous results, we 
have also demonstrated that field theoretical models can yield a neutron skin thickness in 208 Pb 
of less than 0.2 fm. Since a measurement of the neutron radius to about 1% accuracy implies an 
error of about 0.05 fm, the uncertainty in the neutron skin thickness will probably range from 10% 
to 30%, depending upon the magnitude of SR. On this basis, it is expected that the pressure at 
5/8 no could be determined to 20-50%. 

An independent method of constraining S s /S v is from a least squares fit of nuclear models to 
nuclear binding energies. However, this fitting can only reliably establish a correlation between 
S s /S v and S v . We have examined this in the liquid droplet approach using two plausible models. 
The surface tension was parameterized either in terms of (J, n , implying a neutron skin, or in terms 
of 

~ /J-p which allowed for the Coulomb interaction of protons. While these models gave fits to 
the data of similar accuracy, the valley in % 2 differed in the two cases. The behavior of the 
model was shown to be more consistent with the systematic correlation derived from parameter 
optimization of both potential and field-theoretical models of laboratory nuclei. 

The density dependence of the symmetry energy has a number of astrophysical consequences. 
One of the most important is its role in determining the composition of matter. In beta equilibrium, 
the proton fraction is proportional to E sym (n) 3 for small proton fractions. The charge fraction of 
matter plays an essential role in establishing the threshold densities of hyperons and the quark- 
hadron phase transition. If any of these threshold densities are exceeded in a neutron star, the 
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possibility of enhanced neutrino emission by a direct Urea process exists. The threshold for the 
direct Urea process involving nucleons alone is also determined by the charge fraction (for example, 
it is 1/9 in the absence of muons). We have shown that the determination of 5R to have a value 
greater than approximately 0.2 fm would imply that at least the nucleon direct Urea process exists 
in 1.4 M.q neutron stars. 

An additional astrophysical application is the relation we established between 5R and the radii R 
of neutron stars. Although Lattimer and Prakash established a phenomenological relation between 
R for a given mass star and the pressure of neutron star matter above, but near, the nuclear 
saturation density, there exists additional uncertainty in relating this algebraically to the 5R — P 
correlation. However, as shown in this work, a useful relation directly relating 5R and, for example, 
the radius of a 1.4M Q star, R\a, can be established. 

The fact that the neutron skin thickness, which measures symmetry properties below nuclear 
density, and the astrophysical quantities involving the neutron star radius and the Urea threshold 
density, can be correlated is a consequence of a generic trend in the overall density dependence of 
the symmetry energy. This trend was specifically demonstrated by displaying the linear relationship 
which exists between the density derivatives of the symmetry energy at 5no/8 and 3no/2. This 
relationship is independent of the parameterization of the nuclear force. 

Observation of neutron star radii may provide qualitative information on the stellar composi- 
tion. Using the parametric freedom of the EOS's considered, while satisfying the aforementioned 
constraints, we concluded that the minimum radius achievable for a star composed of just nucleons 
and leptons is about 9 km. Observation of a significantly smaller radius would likely imply that 
some softening component (bose condensates, hyperons, quarks, etc.) is present in dense matter, 
although it could also indicate that the EOS's given by potential and field-theoretical models is 
qualitatively incorrect. 

We have also explored the influence of the density dependence of the nuclear symmetry energy 
on the transition pressure marking the boundary between the core and the crust in a neutron star. 
This transition pressure is closely approximated by the phase boundary separating uniform matter 
from matter in which nuclei exist. It was found that the transition pressure was not correlated with 
the neutron skin thickness, although the transition density did show some rather weak correlation. 
More interesting is the fraction of the star's moment of inertia residing in the crust which is 
measurable and which depends on the transition pressure, mass and radius. Here there was some 
correlation with larger values of this fraction tending to imply larger values of the neutron skin 
thickness, but the correlation was not very robust. Probably the neutron skin and the transition 
pressure should be viewed as independent quantities whose measurement would provide separate 
constraints on the EOS. 
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APPENDIX A: MODEL PROPERTIES 



TABLE V: Saturation and surface properties for the models selected in this work. Symbols are uq: nuclear 
matter saturation density in fm~ 3 , B: the binding energy per particle in MeV, K: compression modulus 
in MeV, S v : symmetry energy in MeV, S' v : density derivative of the symmetry energy in MeV fm 3 , <jq: 
surface tension of symmetric matter in MeV fm -2 , 17,5: asymmetry parameter of the surface tension in MeV 
fm~ 2 , and 5R: skin thickness of 208 Pb in fm. The symbol "TW" in the second column refers to new models 
constructed in this work. 
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APPENDIX B: COUPLING STRENGTHS OF MODELS SR2, ES25, ES30, AND ES35 

TABLE VI: Coupling strengths for the model SR2 with nuclear matter equilibrium density tlq = 0.15 fm - , 
binding energy B = 16.38 McV, compression modulus K — 224.5 MeV, Dirac effective mass M*(uq) = 
0.78M, and symmetry energy E sym = 30.1 MeV. The dimensions of a, and bi are such that the Lagrangian 
in Eq. JT3J) is in MeV 4 . 
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11.217 14.055 MeV 0.010091 
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Q2 03 04 
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0.22299 2.2908 xl0~ 7 3.0076 xl0~ 13 


TABLE VII: Coupling strengths for the model es25 with nuclear matter equilibrium density hq = 0.16 fm , 


binding energy 


B = 


16 MeV, compression modulus K 


= 211.7 MeV, Dirac effective mass M*(n ) = 0.76M, 


and symmetry energy E sym = 25 MeV. For the couplings a,; and bi, only the non-zero entries are given. 
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TABLE VIII: Coupling strengths for the model es30 with nuclear matter equilibrium density no — 0.16 fm 3 , 


binding energy 


B = 


16 MeV, compression modulus K 


= 215.4 MeV, Dirac effective mass M*(n ) = 0.66M, 


and symmetry 


enerj 


iy E sym = 30 MeV. Only non-zero 


at and bi are given. 
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0.22445 


TABLE IX: Coupling strengths for the model cs35 with nuclear matter equilibrium density no = 0.1597 fm 3 , 


binding energy 


B = 


16 MeV, compression modulus K 


= 210.0 McV, Dirac effective mass M*(n ) = 0.64AJ, 


and symmetry 


enerj 


IY E sym = 34.93 McV. Only non-zero ai and hi are given. 
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